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airport.  Velocity  and  range  aliasing  will  degrade  the  quality  of  these  estimates,  as 
will  contamination  by  ground  clutter.  For  radars  which  operate  at  short  ranges  and 
at  low  elevation  angles,  as  the  TDWR  will  to  detect  windshears  at  the  airport  sur¬ 
face,  clutter  contamination  is  an  especially  severe  problem.  Multirate  pulse  trains  - 
pulse  trains  containing  multiple  intersample  spacings  -  can  extend  both  the  unambi¬ 
guous  velocity  and  range  of  a  Pulsed  Doppler  Radar  beyond  those  afforded  by  pulse 
trains  with  a  constant  intersample  spacing;  but  the  usual  properties  of  conventional 
clutter  filter  architectures  change  radically  when  applied  to  data  collected  with  a 
multirate  sampling  scheme. 

A  brief  introduction  to  the  systems  and  weather  considerations  for  Doppler 
Weather  Radars  is  provided  and  the  Pulse-Pair  spectral  moment  estimators  are 
presented.  This  introduction  is  followed  by  a  discussion  of  frequency  domain  clutter 
rejection  techniques  for  Batch  PRT  (Pulse  Repetition  Time)  sequences  -  blocks  of 
equispaced  samples  with  the  PRT  alternating  from  block  to  block.  The  mam  topic  of 
the  report  is  clutter  suppression  for  Staggered  PRT  sequences  in  which  the  PRT 
alternates  from  pulse  to  pulse.  The  Staggered  PRT  scheme  has  the  advantage  over 
the  Batch  PRT  scheme  of  spatial  coherency  for  estimates  of  the  radar  return  signal's 
autocorrelation  function  at  the  lags  corresponding  to  the  two  PRT’s.  A  time-varying 
filter  architecture  with  multiple  transfer  functions  is  presented  and  analyzed;  and  its 
interaction  with  the  Pulse-Pair  estimators  is  explored.  Three  design  techniques  for 
Staggered  PRT  filters  are  described  and  assessed  in  the  context  of  clutter  suppres¬ 
sion. 

The  final  section  of  the  report  summarizes  the  results  for  the  Batch  and  Stag¬ 
gered  PRT  schemes  and  provides  suggestions  for  further  research. 
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1.0  INTRODUCTION 

Reliable  weather  parameter  estimates  are  required  of  radars  such  as  the  Ter¬ 
minal  Doppler  Weather  Radar  (TPWR)  -  a  Federal  Aviation  Administration  project 
-  which  will  automatically  detect  hazardous  weather  phenomena  in  the  vicinity  of  an 
airport.  Velocity  and  range  aliasing  will  degrade  the  quality  of  these  estimates,  as 
will  contamination  by  ground  clutter.  For  radars  which  operate  at  short  ranges  and 
at  low  elevation  angles,  as  the  TDWR  will  to  detect  windshears  at  the  airport  sur¬ 
face,  clutter  contamination  is  an  especially  severe  problem.  Multirate  pulse  trains  - 
pulse  trains  containing  multiple  intersample  spacings  -  can  extend  both  the  unambi¬ 
guous  velocity  and  range  of  a  Pulsed  Doppler  Radar  beyond  those  afforded  by  pulse 
trains  with  a  constant  intersample  spacing;  but  the  usual  properties  of  conventional 
clutter  filter  architectures  change  radically  when  applied  to  data  collected  with  a 
multirate  sampling  scheme.  The  topics  of  this  report  are  clutter  suppression  tech¬ 
niques  for  multirate  pulse  sequences  and  their  interaction  with  the  estimation  of 
weather  parameters. 

1.1  BACKGROUND 

The  spatial  scale  for  the  discrete  weather  phenomena  of  interest  in  the  TDWR 
system  ranges  from  tens  of  kilometers,  in  the  case  of  gust  fronts,  to  sub-kilometer 
scales  for  wind  shear  events  such  as  microbursts  or  for  tornadoes.  Sub-kilometer 
resolution  within  50  km  of  the  radar  requires  a  beamwidth  on  the  order  of  1°.  A  1° 
conical  beam  corresponds  to  a  parabolic,  axially-fed,  antenna  whose  diameter  is  10 
m  at  S-Band  i3  GHz)  and  5  m  at  C-band  (6  GHz).  The  requirements  for  angular 
resolution  of  the  weather  radar  serve  to  bound  the  radar  frequency  from  below.  As 
the  radar  frequency  is  decreased  below  S-Band,  it  becomes  increasing  difficult  and 
costly  to  manufacture  both  the  reflector  and  a  pedestal  capable  of  executing 
appropriate  scanning  strategies,  due  to  the  increased  size  and  inertia  of  the 
antenna.  Frequency  bands  below  S-Band  are  not  practical  for  most  ground-based 
weather  radar  applications. 

The  radar  frequency  is  bounded  from  above  by  considerations  of  atmospheric 
attenuation  of  the  radar  signal,  as  it  would  be  unacceptable  for  storms  near  the 
radar  to  block  the  view  of  storms  farther  away  from  the  radar.  If  the  two-way  path 
loss  due  to  atmospheric  attenuation  is  to  be  held  below,  say,  5  dB,  the  radar  fre¬ 
quency  must  not  exceed  7  GHz  [Nathanson69a]  .  So,  weather  radars  for  applica¬ 
tions  such  as  TDWR  are  either  C-Band  or  S-Band  devices.  An  additional  constraint 
on  the  radar  frequency  is  the  availability  of  frequency  allocations  -  C-Band  is 
currently  the  preferred  option  for  TDWR  in  this  context. 

For  Doppler  weather  radars  transmitting  an  equispaced  pulse  train,  the  Pulse 
Repetition  Time  (PRT),  or  intersample  period,  determines  the  radar’s  maximum 
unambiguous  range  ru  and  Nyquist  velocity  v  u.  Denoting  the  speed  of  light  by  c  , 
the  PRT  by  ,  and  the  radar  wavelength  by  X  ,  the  following  expressions  for  ru  and 


1 


v  u  may  be  easily  derived. 
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and  r„ = 


C 


u  2 

Combining  these  two  expressions  to  eliminate  Xj>  ,  we  obtain 

c  X 

vuru- 


(1.1) 


Equation  1.1  represents  the  fundamental  relationship  between  Nyquist  velocity  and 
unambiguous  range  as  a  function  of  radar  frequency.  Figure  1.1  provides  a  graph 
of  this  relationship  for  radar  frequencies  of  3  GHz  and  6  GHz. 


1.2  DOPPLER  POWER  SPECTRUM  AND  PULSE-PAIR  ESTIMATORS 


The  Doppler  power  spectrum  of  a  radar  resolution  cell  represents  the  aggregate 
statistics  of  the  literally  millions  of  hydrom eteors  or  other  particulates  within  that 
volume.  The  average  horizontal  velocity  of  these  particulates  is  approximately  the 
horizontal  windspeed.  A  combination  of  spectral  broadening  mechanisms  including 
shear  within  the  resolution  volume,  turbulence,  scan  modulation  by  the  antenna  pat¬ 
tern  and  the  distribution  of  particle  sizes  results  in  a  Doppler  spectrum  with  Gaus¬ 
sian  shape 
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and  its  Fourier  Transform  partner,  the  weather  autocorrelation  function 
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(1.3a) 


evaluated  at  lag  m ,  where  v  is  the  mean  Doppler  velocity,  av  is  the  spectrum’s 
standard  deviation  or  spectrum  width,  and  S  is  the  power  of  the  weather  echo  pro¬ 
cess. 

The  three  spectral  parameters  of  interest  in  weather  radar  applications  are  the 
first  three  moments  of  the  Doppler  power  spectrum:  the  power  S  ,  which  is  related  to 
the  liquid  water  density  in  the  atmosphere;  the  mean  velocity  v  ;  and  the  spectrum 
width  av  which  is  an  indicator  of  atmospheric  turbulence.  Frequency  domain  esti¬ 
mation  of  these  parameters  is  rarely  employed,  since  the  small  number  of  samples 
per  azimuth  (typically  ~50  )  offers  poor  frequency  resolution  and  since  peak-finding 
algorithms  are  generally  sensitive  to  estimation  errors  resulting  from  measurement 
noise  (which  is  not  included  in  equation  1.2).  Instead,  time  domain  estimation  of  the 
spectral  parameters  is  employed  via  the  autocorrelation  function.  These  estimators 
were  initially  developed  in  the  late  1960’s  [Miller70a]  and  are  commonly  referred  to 
as  pulse-pair  estimators.  One  version  of  these  estimators  makes  use  of  the  auto¬ 
correlation  at  lags  zero  and  one.  only,  and  is  referred  to  as  the  0/1  estimator.  It  is 
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given  by  equations  1.4  [Doviak84aj  . 
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In  these  equations,  Rt  (  0  )  and  Rt(  by)  represent  the  estimated  autocorrelation  at 
lags  zero  and  one  of  the  total  received  signal  -  consisting  of  weather  signal  plus  sys¬ 
tem  noise  -  and  Ns  /  2  is  the  variance  of  the  system  noise.  The  system  noise  is 
assumed  to  be  white  and  independent  of  the  weather  returns.  Ns  /  2  is  a  measur¬ 
able  quantity.  These  estimators  are  popular  for  weather  radar  applications,  as  their 
implementation  is  straightforward. 

Figures  1.2a  and  1.2b  [Doviak84a]  display  the  standard  deviations  of  the  esti¬ 
mators  of  equations  1.4c  and  1.4d,  respectively,  in  the  absence  of  Nyquist  aliasing. 
The  abscissa  is  normalized  spectrum  width,  ov  /  2v  u  ,  the  ordinate  is  the  standard 
deviation  of  the  estimate  normalized  by  the  Nyquist  interval  and  the  number  of 
pulse  pairs  per  dwell  time,  M.  The  curves  in  the  figures  are  parameterized  by 
Signal-to-Noise  ratio  (SNR)  -  10  dB  is  a  representative  figure  for  a  radar  such  as 
the  TDWR. 


As  an  S-Band  example,  consider  the  case  where  \=10  cm  ,  av  =6  m/s,  <^=2  ms 
(av  /2ou  ~  0.25),  an  azimuth  dwell  time  of  50  ms  and  an  SNR  of  20  dB.  Figure 
1.2a  indicates  a  standard  deviation  of  2.0  m/s  -  almost  one-tenth  of  the  Nyquist 
interval  -  for  estimates  of  v  ;  while  Figure  1.2b  indicates  a  standard  deviation  of  1.1 
m/s  for  the  estimate  of  av .  Note  that  with  SNR>0  dB,  these  standard  deviations 
increase  rapidly  as  the  spectrum  width  increases  beyond  0.5  v  u.  Op  should  be 
selected  so  that  v  u  >  2a v  for  all  spectrum  widths  at  which  reliable  Doppler  estimates 
are  required. 


1.3  WEATHER  CHARACTERISTICS 

The  radar  signatures  of  meteorological  phenomena  are  diverse.  Stratiform  pre¬ 
cipitation  may  produce  nearly  uniform  reflectivity  and  velocity  fields  with  small  spec¬ 
trum  width  values.  Tornado  signatures  are  often  evident  in  all  three  fields:  a  hook¬ 
shaped  increased  reflectivity  due  to  spiraling  inflow  of  precipitation  and  airborne 
debris;  high  values  of  spectrum  width;  and  a  velocity  vortex  signature  if  the  tornado 
is  within  50  km  of  the  radar.  The  signatures  of  thunderstorm  outflow  events  such 
as  the  microburst  and  gustfront  are  most  prominent  in  the  velocity  and  spectrum 
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Figure  1.2a:  Standard  deviation  of  pulse-pair  velocity  estimate 
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width  fields  [Uyeda85a]  . 

The  microburst  and  gustfront  are  the  two  most  frequent  types  of  low-altitude 
windshears  at  airports;  their  signatures  must  be  discernible  in  TDWR  Doppler  meas¬ 
urements.  Typical  spectrum  widths  in  microburst  outflows  may  be  on  the  order  of  2 
m/s  [Rinehart87a]  .  The  spectrum  widths  associated  with  gustfronts  are  larger 
than  for  the  microburst  -  in  Appendix  A,  we  recommend  a  design  value  of  6  m/s  in 
this  context.  In  light  of  the  discussion  in  section  1.2,  accurate  (i.e.  standard  devia¬ 
tions  ~  2  m/s)  Doppler  measurements  of  gustfronts  will  require  a  vu  of  at  least  12 
m/s.  This  corresponds  to  minimum  PRF’s  of  480  Hz  and  960  Hz  at  S  and  C-Bands 
(X=0.1  m  and  X=0.05  m),  respectively.  These  minimum  PRF’s  yield  maximum  ru’s 
of  312  km  and  156  km. 

It  has  been  proposed  that  the  TDWR  be  sited  within  15  km  of  its  associated 
airport  and  utilize  elevation  angles  as  low  as  ,'.5°  to  detect  shear  phenomena 
[Merritt87a]  .  Figure  1.3  depicts  the  altitude  coverage  of  a  1°  beam  -  at  an  eleva¬ 
tion  angle  of  0.5°  -  as  a  function  of  ground  range.  Now.  regions  of  significant 
reflectivity  (  e.g.  >  30  dBZ)  in  North  America  may  extend  vertically  to  as  much  as 
14  km  or  45,000  ft  in  altitude  [Vasiloffi34a,  Wilk80a]  .  The  import  of  the  unambi¬ 
guous  range  figure  ru  is  that  the  apparent  range  of  any  target  at  true  range  r,  will 
be  r{modru.  It  is  apparent  from  figure  1.3  that  phenomena  with  significant 
reflectivity  will  be  sensed  beyond  450  km  -  with  a  0.5°  elevation  angle  -  and  that 
they  may  overlay  or  obscure  targets  within  the  unambiguous  range  interval  ru.  The 
former  targets  are  referred  to  as  out-of-trip  targets  or  out-of-trip  weather.  Obscura¬ 
tion  by  out-of-trip  weather  severely  limits  the  radar’s  detection  capability  for  in-trip 
hazards  with  low  reflectivity  such  as  gustfronts  or  microbursts  in  the  American 
high-plains  regions  [Fujita85a]  . 

The  minimum  PRF’s  above,  based  on  considerations  of  estimate  accuracy,  are 
not  low  enough  to  avoid  obscuration  by  out-of-trip  weather.  Nor  are  they  high 
enough  to  ensure  unambiguous  Doppler  velocity  estimates.  Peak  microburst 
windspeeds  of  33  m/s,  for  example,  have  been  observed  [Fujita85a]  ;  and 
windspeeds  associated  with  tomadic  activity  may  be  as  high  as  90  m/s 
[Doviak84a]  .  Pulsed  Doppler  Weather  Radars  operating  at  S  and  C-Bands  must 
contend  with  degrees  of  both  range  obscuration  and  velocity  aliasing. 

1.4  RANGE  DE-ALIASING 

Range  obscuration  is  difficult  to  predict  in  advance  of  actual  measurements; 
and  no  satisfactory  techniques  exist  for  removing  the  effects  of  obscuration  by  post¬ 
processing  of  the  measured  data.  Adaptive  PRT  selection  and  random  phase  coding 
of  the  pulse  train  are  two  methods  for  combatting  obscuration. 

In  the  former  case,  a  rapid  surveillance  scan  with  a  large  PRT  (large  ru)  is 
employed  to  unambiguously  acquire  regions  of  reflectivity.  This  information  is  used 
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Figure  1.3:  Vertical  profile  of  1°  beam  at  0.5°  elevation  versus  ground  range 
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to  select  the  PRT  for  a  subsequent  weather  collection  scan.  Random  phase  coding, 
on  the  other  hand,  whitens  the  received  signal  component  due  to  the  out-of-trip 
weather,  at  the  expense  of  decreasing  the  effective  SNR  for  the  in-trip  weather. 

1.5  VELOCITY  DE-ALIASING 

The  need  to  free  radar  velocity  data  from  ambiguities  is  especially  critical  if  the 
data  are  to  be  processed  by  automated  algorithms.  Unlike  range  obs deration,  how¬ 
ever.  techniques  do  exist  for  "de-aliasing"  aliased  velocity  estimates. 

There  are  two  general  approaches  to  velocity  de-aliasing.  The  first  approach, 
relative  de-aliasing,  invokes  assumptions  of  local  radial  and  azimuthal  continuity  in 
processing  the  measured  velocity  field.  As  an  example,  imagine  a  radar  with 
uu  =  10  m/s.  Suppose  that  a  radial  of  velocity  measurements  has  just  become 
available.  Suppose  further  that  the  first  hundred  range  gates  have  velocity  values 
of  S  to  10  m/s,  and  that  the  101st  gate’s  velocity  value  is  -9  m/s.  The  logic  of 
radial  continuity  would  declare  that  a  velocity  alias  had  occurred  and  reset  the  value 
of  -9  m/s  to  11  m/s.  Azimuthal  continuity  operates  in  an  analogous  fashion. 

Relative  de-aliasing  produces  a  velocity  field  containing  extended  regions  which 
differ  from  the  "true"  Doppler  velocity  field  by  an  additive  constant  equal  to  some 
multiple  of  the  Nyquist  interval  length,  2vu.  The  second  approach  to  velocity  de¬ 
aliasing,  absolute  de-aliasing,  attempts  to  generate  an  estimate  of  the  "true" 
Doppler  windfield.  This  can  be  achieved  by  the  use  of  pulse  trains  with  multiple 
PPT’s.  The  mechanism  of  absolute  de-aliasing  will  be  discussed  subsequently. 

For  the  identification  of  phenomena  whose  discriminant  is  spatial  variations  in 
the  velocity  field,  relative  de-aliasing  of  that  field  is  sufficient;  relative  de-aliasing 
produces  a  datafield  whose  spatial  "gradient"  locally  matches  that  of  the  true 
Doppler  velocity  field.  This  class  of  phenomena  includes  the  microburst  and  other 
shear  events.  Absolutely  de-aliased  velocity  estimates  are  required  when  the  velo¬ 
city  data  is  used  directly  -  as  in  the  Sectorized  Uniform  Wind  Algorithm  [Witt87a]  , 
which  is  used  to  estimate  the  windshift  associated  with  the  passage  of  a  gustfront. 

The  discussions  in  the  remainder  of  this  report  will  concentrate  on  issues  related 
to  absolute  de-aliasing. 

1.6  ABSOLUTE  VELOCITY  DE- ALIASING 

Section  1.2  described  the  estimation  of  a  discrete  mean  Doppler  velocity  when 
the  radar  pulse  train  contains  a  single  PRT.  Now,  suppose  that  the  weather  process 
is  wide-sense  stationary  over  an  azimuth  dwell-time  -  a  reasonable  assumption  so 
long  as  the  resolution  volume  contains  a  homogeneous  windfield  -  and  that  we  are 
able  to  independently  interrogate  the  weather  process  with  two  different  PRT’s,  c  1<Sr 
and  c2& r,  where  cx  and  c2  are  relatively  prime,  positive,  integers  such  that  Cj<c2. 
The  Nyquist  velocities  corresponding  to  these  two  sets  of  measurements  are 
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i>ul  =  X(4c1Xr)~1  and  vu2  =  X (  4c2<t )_1.  respectively.  Generating  separate  mean 
velocity  estimates  (  e.g.  by  using  1.4c  twice)  we  obtain  two.  potentially  aliased,  esti¬ 
mates  v  x  and  v  2. 

and  v2  may  be  combined  to  generate  a  de-aliased  estimate  with  the  higher 
effective  Nyquist  velocity  of  v  u  ,  the  Nyquist  velocity  corresponding  to  a  PRT  of  Xf. 
This  statement  may  seem  counterintuitive,  but  it  derives  directly  from  number 
theory  which  tells  us  that  the  ordered  pair  {  (u  mod  vu])  ,  ( v  mod  vu2 )  }  is  a  unique 
function  of  v  up  to  the  Least  Common  Integral  Multiple  of  uu  x  and  uu2  which  is 
civui  =c2vu2  =  vu  =  M4&j<)-1  -  An  example  of  this  is  provided  in  Figure  1.4  for 
the  case  c  x  =2  ,  c2  =3.  The  top  pair  of  curves  in  the  figure  display  the  Doppler  fre¬ 
quency  abasing  for  the  two  sampling  frequencies  F1=(c]Xr)~1  Hz  and 
F2  =  (c2^r)_1  Hz.  The  bottom  curve  displays  the  difference  of  the  two  abased  fre¬ 
quency  estimates,  f  5  and  f  2  as  a  function  of  true  Doppler  frequency.  Note  that  this 
difference  is  periodic  with  period  Xf 1  ,  the  Nyquist  interval  corresponding  to  a  PRT  of 
Xr  • 

Thus,  if  we  can  "simultaneously"  interrogate  the  weather  process  with  two 
PRT’s  c iXr  and  c2 Xp  ,  we  can  obtain  an  estimate  with  a  Nyquist  velocity  of 
X  ( 4XJ-)-1  by  comparing  the  estimates  v  x  and  v2,  and  determining  the  appropriate 
additive  correction  (  e.g.  via  a  table  look-up).  This  final  step  in  unfolding  v  x  and  v2 
is  a  statistical  detection  problem  which  is  sensitive  to  errors  in  those  two  estimates. 
For  the  example  of  Figure  1.4,  errors  of  magnitude  greater  than  'AF  2  in  f  x—  f2 
would  result  in  an  erroneous  unfolded  estimate  (assuming  that  the  de-aliasing  algo¬ 
rithm  examines  only  the  difference  f  x—f  2  ). 

Pulse  trains  which  incorporate  two  PRT’s  are  of  two  basic  types.  The  first  is  a 
Batch  PRT  sequence  in  which  half  the  azimuth  dwell  time  is  spent  taking  samples 
at  PRT  CjXp  and  the  other  half  is  spent  taking  samples  at  PRT  c2Xp  .  The  second 
type  of  pulse  train  is  a  Staggered  PRT  sequence  in  which  the  PRT  alternates  from 
CjXj-  to  c2Xp  on  a  pulse-by-pulse  basis.  Both  types  of  sequences  will  be  discussed  in 
more  detail.  We  note,  though,  that  the  Batch  PRT  scheme  is  prone  to  difficulties 
when  shears  exist  within  the  resolution  volume  -  a  likelihood  at  longer  ranges  -  for 
the  following  reason.  The  logic  of  velocity  de-abasing  assumes  that  v  x  and  u  2  are 
differently  abased  estimates  of  the  same  mean  velocity.'  If  the  windfield  in  the  reso¬ 
lution  volume  is  not  uniform,  the  v  }  and  v2  generated  with  a  Batch  PRT  scheme  wifi 
be  abased  estimates  of  the  mean  velocities  for  statistically  different  weather 
processes.  In  this  case,  the  assumptions  of  velocity  de-abasing  will  be  violated  and 
the  de-abasing  technique  will  generally  fail.  De-abasing  failures  are  not  graceful 
since  the  errors  are  on  the  order  of  the  Nyquist  intervals  corresponding  to  the  two 
PRT’s.  The  Staggered  PRT  scheme  is  a  better  choice  in  the  respect  that  v  x  and  v2 
will  be  spatially  coherent.  Some  results  on  the  performance  of  Staggered  PRT  velo¬ 
city  estimation  are  provided  in  [Sirmans76a]  . 
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VELOCITY  DE-ALIASING  WITH  MULTIPLE 

SAMPLE  RATES 


TWO  PRFs:  F.  =  .  F,  =  -1.  ,  dT  IS  "BASE  INTERVAL" 

CA  C2dT 


WHEN  C1  AND  C2  ARE  RELATIVELY  PRIME  INTEGERS.  THE  EFFECTIVE 

NYQUIST  FREQUENCY  IS  F.  =  -i-  _ 

Z0T  e  g.  C1  =  2.  C2  =  3 


Figure  1.4:  Velocity  de-aliasing  with  multiple  sample  rates 


11 


1.7  THE  CLUTTER  CHALLENGE 


The  final  topic  to  be  discussed  in  this  introduciion  is  the  need  for  effective  clutter 
suppression  in  Doppler  weather  radars.  For  land  based  weather  radars,  clutter  con¬ 
sists  primarily  of  returns  from  ground  targets  such  as  trees,  buildings  and  oro¬ 
graphic  features.  Unfortunately,  these  targets  possess  radar  cross-sections  which 
are  typically  much  larger  than  the  aggregate  radar  cross-section  of  the  particulates 
within  a  radar  resolution  volume. 

Clutter  of  a  given  radar  cross-section  assumes  different  equivalent  dBZ  values 
as  a  function  of  range.  Recent  results  [Mann86a]  have  shown  that  stressful  clutter 
environments  correspond  to  effective  clutter  returns  of  55  dBZ  at  10  km  range  to  45 
dBZ  at  100  km  range,  while  average  clutter  returns  are  on  the  order  of  45  dBZ  at 
10  km  range  to  35  dBZ  at  100  km  range.  If  we  wish  to  maintain  a  10  dB  signal  to 
clutter  ratio  for  low- reflectivity  phenomena  such  as  gustfronts.  clutter  suppression  in 
excess  of  50  dB  is  required.  With  antenna  rotation  rates  on  the  order  of  4  rpm,  this 
suppression  must  be  achieved  within  azimuth  dwell  times  on  the  order  of  40  ms.  We 
note  that  the  azimuth  dwell  time  is  an  important  system  parameter  for  clutter 
suppression.  Its  inverse  is  proportional  to  the  attainable  transition  bandwidth  of 
FIR  clutter  filters  that  operate  on  equispaced  data,  and  to  the  frequency  resolution 
available  for  frequency-domain  clutter  suppression  techniques. 

1.8  SUMMARY  OF  THE  ISSUES  AND  OUTLINE  OF  THE  REPORT 

We  have  explained  why  Doppler  weather  radars  usually  operate  in  the  S-band 
or  C-Band  frequency  regmies  (2-7  GHz),  and  that  accurate  Doppler  estimates  in  the 
context  of  TDWR  necessitate  a  minimum  Nyquist  velocity  of  12  m/s  -  PRF's  of  480 
Hz  and  960  Hz  at  X=10  cm  and  \=5  cm,  respectively.  We  have  shown  that  these 
PRF’s  are  neither  low  enough  to  prevent  range  obscuration,  nor  high  enough  to 
prevent  velocity  aliasing.  Schemes  for  de-aliasing  velocity  estimates  were  described; 
and  we  noted  that  it  is  possible  to  absolutely  de-alias  velocity  estimates  by  use  of 
two  PRT’s.  Two  types  of  radar  pulse  trains  incorporating  multiple  PRT’s  were 
described,  the  Batch  PRT  and  Staggered  PRT  waveforms.  We  have  also  showm  that 
effective  clutter  suppression,  in  excess  of  50  dB,  is  required  if  shear  events  with  low 
reflectivities,  such  as  gustfronts,  are  to  be  detected. 

The  organization  of  the  remainder  of  this  report  is  outlined  below. 


§2  Clutter  Suppression  for  Batch  PRT  Sequences 

The  rationale  for  performing  Batch  PRT  clutter  suppression  in  the  frequency 
domain  is  laid  out  and  moment  estimation  wdth  this  scheme  is  described. 

§2.1  Data  Taper  Selection  for  Frequency  Domain  Clutter  Suppression 

A  normalized  expression  for  assessing  data  taper  performance  in  the  con¬ 
text  of  clutter  suppression  is  derived.  This  expression  is  employed  to 
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evaluate  seven  candidate  windows  for  Batch  and  equispaced  PRT  clutter 
suppression  in  a  TDWR  context. 

§3  Clutter  Suppression  for  Staggered  PRT  Sequences 
§3.1  Staggered  PRT  Filter  Structure 

The  staggered  PRT  sampling  scheme  is  described  and  the  architecture  of 
the  Staggered  PRT  Filter  is  introduced. 

§3.2  Staggered  PRT  Filter  Response 

Z-transforms  for  multirat*  data  sequences  are  presented  and  used  to  jus¬ 
tify  a  time- varying  filter  architecture.  The  effects  of  having  identical  par¬ 
tial  response  sequences,  which  is  equivalent  to  a  conventional  FIR  tapped 
delay  line,  are  discussed  ff  :*  odd  and  even  filter  lengths. 

§3.2.1  Average  Pow'er  Transfer  Function 

The  staggered  PRT  filter’s  power  transfer  function  for  lags  which  are 
multiples  of  cx+c2'  |G(~)  |".  is  derived.  This  is  one  of  the  three  stag¬ 
gered  PRT  filter  transfer  functions  that  interact  with  the  pulse-pair 
moment  estimators.  Multiple  transfer  functions  are  an  important 
feature  of  the  time-varying  filter  architecture. 

§3.2.2  Convolutions  of  Staggered  PRT  Sequences 

The  convolution  of  staggered  PRT  sequences  is  contrasted  with  the  con¬ 
volution  of  equispaced  sequences  as  preparation  for  section  3.2.3. 

§3.2.3  Lag  and  Lag  c2  Power  Transfer  Functions 

The  staggered  PRT  filter's  power  transfer  functions  for  lags  of  c,  or  c2 
plus  multiples  of  c  i  -tc  2  are  derived. 


§3.3  Properties  of  the  C,  (z  )  and  Effective  Power  Transfer  Functions 
§3.3.1  Z-Plane  Zero  Distribution  of  |G(u) 


Conjugate  and  reciprocal  symmetry  of  the  z-plane  zero  distribution  of 
|G(-v)  | 2  is  established.  These  symmetries  guarantee  a  IGM  |2  which 
is  a  real,  even,  function  of  frequency. 


§3.3.2  Phase  Responses  of  C  (*  (e',u  T)C  and  C 0(e;  ‘  T1C  *  (eJ'“‘T) 

The  phase  response  of  the  lags  Cj  and  c2  power  transfer  functions  are 
discussed.  The  non-zero  phase  response  of  the  staggered  PRT  filter  is  a 
key  factor  in  its  interaction  with  the  moment  estimation  process. 
Examples  are  provided  of  velocity  estimation  with  a  staggered  PRT 
filter. 


§3.3.3  Uncontrollable  Zeroes  in  the  C,(z  ) 

Not  all  of  the  z-plane  zero  locations  of  the  staggered  PRT  filter  may  be 
directly  specified  by  the  filter  designer.  The  number  and  z-plane  distri¬ 
bution  of  these  uncontrollable  zeroes  are  discussed  as  functions  of  the 
filter  length  and  stagger  ratio. 
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^4  Staggered  PRT  Filter  Design  Techniques 


^4.1  Design  Technique:  Coefficient  Selection 

A  design  technique  bused  on  selecting  coefficients  from  a  time-aligned 
equispaced  filter  with  a  PRT  of  I'Vp  is  described  and  assessed.  The  DFT  is 
shown  to  be  an  isometric  mapping  between  the  time  and  frequency 
domains. 

$4.2  Design  Technique:  Maximally  Flat  Response 

The  assignment  of  all  controllable  zeroes  to  z=l  is  assessed  as  a  design 
technique. 

$4.3  Design  Technique:  Output  Signal  to  Clutter  Ratio  Maximization 

Staggered  PRT  filter  coefficients  are  selected  to  provide  optimal  suppres¬ 
sion  of  scan- modulated  clutter  with  a  generalized  weather  spectrum. 

$5  Conclusions 

ml  Batch  PRT 
$.~>.2  Staggered  1  MIT 
S~).3  Further  Research 

Suggestions  are  provided  for  characterizing  the  z-plane  zero  locations  of 
the  staggered  PRT  filter;  and  techniques  for  phase  control  of  the  lags  C) 
and  c_.  transfer  functions  are  described.  A  composite  filter  architecture, 
combining  equispaced  and  staggered  PRT  filters,  is  described. 
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2.0  CLUTTER  SUPPRESSION  FOR  BATCH 
PRT  SEQUENCES 

The  Batch  PRT  pulse  sequence  is  a  pulse  train  with  two  alternating  blocks  of 
equispaced  pulses  as  depicted  in  Figure  2.1.  A  block  of  N  x  samples  is  collected  with 
PRT  cxt>r  ,  followed  by  a  block  of  N2  samples  with  PRT  c2<^r  •  Nx  and  N2  are  typi¬ 
cally  selected  so  that  each  block  of  equispaced  samples  has  a  duration  of  approxi¬ 
mately  half  the  azimuth  dwell  time.  Clutter  rejection  is  performed  independently  for 
each  block  of  equispaced  samples. 

The  Doppler  spectrum  of  the  weather  process  may  be  concentrated  near  zero 
velocity.  As  it  is  undesirable  for  a  clutter  filter  to  reject  the  weather  signal,  the 
clutter  filter’s  stopband  and  transition  regions  should  be  limited  to  the  narrowest 
velocity  interval,  centered  about  zero,  which  affords  the  necessary  rejection.  At  S- 
Band.  typical  spectrum  widths  for  scan-modulated,  distributed  clutter  are  ~  0.25 
m/s  with  the  assumption  of  a  Gaussian  shaped  clutter  spectrum  [Evans83a]  .  For 
this  model,  50  dB  clutter  rejection  requires  that  the  filter’s  stopband  extend  to  at 
least  four  standard  deviations,  a  1  m/s  stopband  edge.  The  transition  region 
bandwidth  should  be  less  than  2  m/s,  so  that  passband  performance  is  assured  for 
velocities  greater  than  3  m/s. 

For  FIR  clutter  filters  operating  upon  equispaced  data  ( e.g .  passband  power 
ripple  of  3  dB,  stopband  power  ripple  below  -50  dB,  and  more  than  10  taps),  the 
width  of  the  transition  region,  in  Hz,  is  approximately  1.35  times  the  inverse  of  the 
filter’s  impulse  response  duration  in  seconds  [Herrmann73a]  .  Typical  antenna 
rotation  rates  for  weather  radars  correspond  to  dwell  times  on  the  order  of  50  ms. 
The  dwell  time  per  equispaced  block  of  samples  is  then  25  ms;  and,  if  we  restrict  the 
filters  impulse  response  duration  to,  say,  15  ms  to  insure  a  sufficient  number  of 
samples  at  the  filter’s  output  for  reliable  estimates  of  the  autocorrelation  at  lags  zero 
and  one,  the  resulting  transition  region  bandwidth  is  90  Hz,  or  4.5  m/s  at  S-Band. 
This  example  provides  an  optimistic  estimate  of  the  transition  region  bandwidth  - 
the  10  ms  integration  time  may  not  provide  pulse-pair  Doppler  moment  estimates  of 
sufficiently  small  variance  [Zmic79a]  . 

Rather  than  performing  clutter  suppression  in  the  time  domain,  weather  radars 
employing  Batch  PRT  sequences  could  perform  clutter  rejection  in  the  frequency 
domain  [Passarelli81a]  .  Each  block  of  equispaced  samples  is  processed  in  the  fol¬ 
lowing  manner.  The  data  sequence,  possibly  weighted  by  a  data  taper,  is 
transformed  into  the  frequency  domain  via  the  Discrete  Fourier  Transform  (DFT). 
The  DFT  of  the  transformed  sequence  is  then  replaced  at  each  point  by  its  magni¬ 
tude  squared,  generating  a  periodogram  estimate  of  the  weather  plus  clutter  spec¬ 
trum.  The  centra!  (zero-frequency)  coefficient  of  the  periodogram  is  set  to  zero,  pos¬ 
sibly  along  with  its  neighbors;  and  then  transformed  back  via  the  Inverse  DFT  to 
generate  an  estimate  of  the  weather  autocorrelation  function.  Pulse-pair  estimation 
is  then  performed  with  the  autocorrelation  estimates. 
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Figure  2.1: 


TIME  (  6Tsec) 


Batch  PRT  sampling  scheme 


2.1  DATA  TAPER  SELECTION  FOR  FREQUENCY  DOMAIN  CLUTTER 
SUPPRESSION 


The  particular  data  taper  which  is  applied  to  time  domain  sequence  prior  to 
performing  the  DFT  has  a  significant  effect  on  the  efficacy  of  frequency-domain 
clutter  suppression.  For  periodogram-based  spectral  estimation,  data  tapers  are 
employed  to  limit  the  blurring  or  apodization  of  the  spectral  estimate  resulting  from 
a  finite  record  length  [Harris78a]  .  In  an  analogous  fashion,  data  tapers  may  be 
employed  for  frequency-domain  clutter  suppression  to  confine  the  clutter’s  spectral 
energy  to  Doppler  frequencies  near  zero. 

The  intent  of  this  current  section  is  to  present  an  analytical  expression, 
independent  of  radar  constants,  which  is  useful  in  assessing  the  effectiveness  of  vari¬ 
ous  data  tapers  for  frequency-domain  clutter  suppression.  Assume  a  finite  energy, 
wide-sense  stationary,  clutter  process  s  (t )  whose  true  power  spectrum  is  given  by 


Ss  (/ ) 


A  *0  „-f2/ 2,jc2 

Tn 


where  oc  is  the  clutter  spectrum  width  in  Hz.  We  will  denote  the  autocorrelation 
function  "f  a  it )  by  Rs  (r).  Suppose  that  we  make  an  observation  of  s  it )  of  duration 
T  and  weight  that  observation  by  a  data  taper 


wit)  = 


d«>  .  \t  I  <  ~ 

° ,  i< i>£ 


where  d  it )  is  a  real,  even,  time  function  with  finite  energy. 

Now,  the  continuous  analogue  of  the  windowed  periodogram  estimate  is  the  fol¬ 
lowing  expectation  with  respect  to  s  (t ). 
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J  wit)s  ( t  )e  d t 
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(2.1) 


Since  w  ( t )  and  s  ( t )  both  possess  finite  energy,  the  integrals  may  be  combined  to 
yield 


w  i^+r)w  *  (t)s  (£+r)s  *  {r)e  y  2^ f  d£dr 


where  we  have  made  the  substitution  £  =  t  —  r  .  Noting  that  w'  {t)=w(—  '),  we  may 
take  the  expectation  inside  the  integral  and  change  the  order  of  integration  with  t 
replaced  by  -r  to  obtain 
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where  W (f  )  is  the  Fourier  Transform  of  the  data  taper  w  ( t ). 

Now  suppose  that  we  set  this  "periodogram"  to  zero  for  frequencies  in  the  inter¬ 
val  [— ft  ,  17J  .  As  W(f  )  and  5S  (f  )  are  both  even,  the  power  remaining  in  the  spec¬ 
trum  is  given  by 

OO  OC 

Ep  ^2  /  J  W2{OSs(f- Od&f  (2.2) 

n  — oo 

OO  00 
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where  we  have  made  the  replacement  q  =  - — — .  Substituting  for  ,  we  have 
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yielding  our  result 

OO 

Ep  =  2K0  J  W2(0  erfc 

~O0 

where  erfc(x )  is  the  complimentary  error  function, 

1  T  -v 

erfc(x )  =  1 - f  e  2  d<,* 

V2tt  _m 

Equation  2.3  is  a  useful  result  in  that  it  completely  characterizes  the  efficacy  of 
a  data  taper  for  frequency-domain  clutter  rejection,  parameterized  by  the  ratio  of  the 
clutter  spectrum  width  ac  and  the  stopband  cut-off  frequency  17  .  We  can  apply  this 
result  to  assess  frequency-domain  clutter  rejection  for  blocks  of  equispaced  pulses  by 
plotting  Ep  as  a  function  of  17/crc  for  a  sampled  data  taper  wk  of  length  N  and 
azimuth  dwell  time  N  . 

Figure  2.2  displays  normalized  graphs  of  equation  2.3  for  seven  different  data 
tapers,  with  representative  S-Band  Doppler  parameters.  This  example  is  for  the 
case  of  a  uniform  sampling  frequency  of  640  Hz  over  the  azimuth  dwell  time  and  a 
33  point  data  taper.  The  azimuth  dwell  time  is  52  ms.  The  value  selected  for  oc  is 
5  Hz,  which  corresponds  to  a  0.25  m/s  clutter  spectrum  width  for  a  3  GHz  radar. 


Rs(0  f 


-r)w(r)e  drdf  = 


i  S,(f)  *  W2(f)  } 
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PASSBAND  EDGE  I  CLUTTER  WIDTH 


SMPL.  FREQ.  =  640.00  Hz 
CLUTTER  WIDTH  =  5.0000  Hz 
33  POINT  WINDOW 
A  =  DIRICHLET  WINDOW 
B  =  BARTLETT  WINDOW 


C  =  PARZEN  WINDOW 

D  =  HAMMING  WINDOW 

E  =  HANN  WINDOW 

F  =  BLACKMAN  WINDOW 

G  =  DOLPH-CHEBYSHEV  WINDOW 


Figure  2.2:  Clutter  rejection  via  the  windowed  periodogram,  33  point  window 
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Window  "A"  is  the  Dirichlet  or  "boxcar"  window  which  applies  uniform  weights 
to  the  data  sequence  -  the  Dirichlet  window  corresponds  to  the  case  in  which  no  data 
taper  is  employed.  Window  "B”  is  the  Bartlett  or  triangular  window,  it  is  the  time- 
scaled  convolution  of  the  Dirichlet  window  with  itself.  Window  "C"  is  the  the  Parzen 
window,  it  is  the  time-scaled  convolution  of  the  Barlett  window  with  itself  and  con¬ 
sists  of  segments  of  second  and  third  order  polynomials.  Window  "D",  the  Hamming 
window,  consists  of  a  cosine  bell  upon  a  pedestal;  it  is  a  weighted  average  of  Window 
"E"  the  von  Hann  window,  which  is  a  simple  cosine  bell,  and  the  Dirichlet  window. 
The  Blackman  window,  "F",  is  a  second  order  trigonometric  polynomial.  The  Dolph- 
Chebyshev  window,  "G",  has  uniform  sidelobe  levels,  and  possesses  the  narrowest 
main  lobe  width  for  a  given  sidelobe  level  (The  examples  in  this  report  employ  a 
squared  transform  sidelobe  level  of  -70  dB,  this  value  was  selected  to  provide  50  dB 
of  clutter  rejection  in  the  example  of  Figure  2.2.)  [Harris78a,  Jasik61a]  . 

These  seven  windows  provide  a  representative  sample  of  data  tapers;  ranging 
from  the  Dirichlet  window  which  possesses  a  very  narrow  main  lobe  in  the  frequency 
domain  and  large  sidelobes,  to  the  Blackman  window  which  possesses  sidelobes 
which  are  orders  of  magnitude  less  than  those  of  the  Dirichlet  window  (at  the 
expense  of  increased  main  lobe  width),  to  the  Chebyshev  window  which  "optimally" 
balances  main  lobe  width  against  sidelobe  level.  The  dynamic  ranges  of  the  win¬ 
dows’  coefficients  are  0,  12,  33,  11,  20,  25  and  18  dB  for  windows  A  through  G  in 
Figure  2.2,  respectively.  Windows  A,  B  and  C  are  "constructed"  windows;  while  win¬ 
dows  D,  E  and  F  are  trigonometric  windows.  Window  G’s  DFT  is  a  sampled  Che¬ 
byshev  polynomial. 

In  examining  Figure  2.2,  the  frequency  bin  size  of  the  DFT  must  be  considered 
(  i.e.  the  frequency  locations  at  which  coefficients  may  be  set  to  zero).  In  this  case 
the  bin  size  is  (640/33)=19.4  Hz  »  1  m/s  at  S-Band.  It  is  apparent  that  some  sort 
of  data  taper  is  required  for  frequency-domain  clutter  rejection,  as  curve  A  -  the 
taperless  case  -  fails  to  reach  even  20  dB  suppression  with  a  passband  edge  of  3 
m/s  (Q/ac  —  12)  .  The  knees  in  the  various  curves  correspond  to  transitions  from 
mainlobe  to  sidelobe  rejection  for  the  windows.  This  transition  happens  at  the  first 
non-zero  DFT  bin  edge  for  window  A  (  Q/ ac  =4  s  20  Hz  or  1  m/s  );  at  the  second 
non-zero  DFT  bin  edge  for  windows  B,  D  and  E;  and  at  the  third  non-zero  DFT  bin 
edge  for  windows  C,  F  and  G.  This  figure  demonstrates  that  acceptable  (  i.e.  >  50 
dB  )  clutter  rejection  is  achievable  with  this  technique,  for  a  3  m/ s  passband  edge  at 
S-Band.  We  may  also  infer  from  the  figure  that  the  maximum  achievable  clutter 
suppression  is  driven  by  a  data  taper’s  sidelobe  levels  rather  than  its  main  lobe 
width. 

The  parameters  for  Figure  2.3  are  identical  to  those  of  Figure  2.2  except  that 
the  azimuth  dwell  time  has  been  reduced  to  40  ms.  Reducing  the  dwell  time  from 
50  to  40  ms  has  the  approximate  effect  of  scaling  the  curves  in  Figure  2.3  horizon¬ 
tally  by  a  factor  of  50/40  relative  to  Figure  2.2.  Since  each  DFT  point  nominally 
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SMPL.  FREQ.  =  640.00  Hz  C  =  PARZEN  WINDOW 
CLUTTER  WIDTH  =  5  0000  Hz  D  =  HAMMING  WINDOW 
25  POINT  WINDOW  E  =  HANN  WINDOW 

A  =  DIRICHLET  WINDOW  F  =  BLACKMAN  WINDOW 

B  =  BARTLETT  WINDOW  G  =  DOLPH-CHEBYSHEV  WINDOW 


Figure  2.3:  Clutter  rejection  via  the  windowed  periodogram,  25  point  window 


corresponds  to  a  frequency  band  of  one  bin  width,  centered  at  that  point,  it  appears 
that  50  dB  of  rejection  with  a  40  ms  dwell  time  could  be  realized  by  zeroing  the 
same  number  of  DFT  coefficients  as  for  a  50  ms  dwell  time. 

Figure  2.4  is  an  S-Band  example  of  frequency-domain  clutter  suppression  with 
a  Block  PRT  scheme.  All  parameters  are  identical  to  those  of  Figure  2.2,  except 
that  the  window  length  has  been  halved,  corresponding  to  a  block  of  equispaced  sam¬ 
ples  which  occupy  half  an  azimuth  dwell  time  as  in  Figure  2.1  .  The  DFT  bin  size 
increases  accordingly  to  (640/ 1 7 )=37 .6  Hz  *  2  m/s  .  The  available  passband  edges 
are  then  at  2,4,6,  ...  m/s.  By  halving  the  observation  time,  the  spread  of  the  clutter 
spectral  energy  due  to  a  finite  record  length  is  increased,  roughly  by  a  factor  of  two, 
and  clutter  rejection  for  a  given  passband  edge  is  correspondingly  reduced.  Figure 
2.2  displays  average  rejection  of  24  dB  for  a  passband  edge  of  2  m/s  (  Q/ ac  =8  ) , 
while  the  corresponding  rejection  depicted  in  Figure  2.4  is  closer  to  10  dB. 

Figure  2.5  displays  the  available  clutter  rejection  for  the  case  of  Figure  2.4,  but 
for  passband  edges  up  to  8  m/s  (  Q/ac  =  32  ).  Figure  2.5  closely  resembles  Figure 
2.2  with  the  abscissa  scaled  by  a  factor  of  two.  The  requisite  50  dB  of  rejection  is 
attained  with  the  Blackman  and  Dolph-Chebyshev  data  tapers,  but  not  until  the 
passband  edge  is  moved  to  the  third  DFT  bin  at  6  m/s.  Relaxing  the  passband  edge 
requirement  from  3  m/s  to  6  m/s  will  not  necessarily  make  this  an  acceptable 
clutter  rejection  scheme  for  use  with  velocity  de-aliasing.  With  a  6  m/s  passband 
edge,  the  filter  stopband  occupies  three-eighths  of  the  entire  Nyquist  interval  for  this 
case.  Not  only  will  weather  returns  with  Doppler  velocities  in  the  range  [-6  m/s  ,  6 
m/s]  be  suppressed;  velocities  in  the  Doppler  intervals  which  alias  to  that  range, 
[  —  6+2nu  u  ,  6+2m;  u  ]  where  n  el,  will  also  be  suppressed. 

If  frequency-domain  clutter  suppression  is  to  be  employed  with  Batch  PRT 
schemes  for  velocity  de-aliasing  purposes,  one  of  the  following  two  compromises  must 
be  made.  Either 

1)  the  radar’s  azimuth  dwell  time  must  be  increased  beyond  50  ms  - 
equivalently,  the  azimuth  scanning  rate  must  be  reduced  below  20 
degrees  per  second  (increasing  the  likelihood  of  velocity  de-aliasing 
errors  due  to  temporal/  azimuthal  non-stationarity  of  the  weather 
process);  or 

2)  the  50  dB  clutter  rejection  requirement  must  be  relaxed,  at  the  cost 
of  increased  potential  for  obscuration  of  low  reflectivity,  low-altitude, 
phenomena. 
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NORMALIZED  POST-FILTER  POWER  (dB) 


PASSBAND  EDGE  /  CLUTTER  WIDTH 

SMPL.  FREQ.  =  640.00  Hz  C  =  PARZEN  WINDOW 
CLUTTER  WIDTH  =  5.0000  Hz  D  =  HAMMING  WINDOW 
17  POINT  WINDOW  E  =  HANN  WINDOW 

A  =  DIRICHLET  WINDOW  F  =  BLACKMAN  WINDOW 

B  =  BARTLETT  WINDOW  G  =  DOLPH-CHEBYSHEV  WINDOW 


Figure  2.4:  Clutter  rejection  via  the  windowed  periodogram,  17  point  window 


PASSBAND  EDGE  /  CLUTTER  WIDTH 


SMPL.  FREQ  =  640.00  Hz 
CLUTTER  WIDTH  =  5.0000  Hz 
17  POINT  WINDOW 
A  =  DIRICHLET  WINDOW 
B  =  BARTLETT  WINDOW 


C  =  PARZEN  WINDOW 

D  =  HAMMING  WINDOW 

E  =  HANN  WINDOW 

F  =  BLACKMAN  WINDOW 

G  =  DOLPH-CHEBYSHEV  WINDOW 


Figure  2.5:  Clutter  rejection  via  the  windowed  periodogram,  17  point  window 
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3.0  CLUTTER  SUPPRESSION  FOR  STAGGERED 

PRT  SEQUENCES 

Clutter  suppression  for  staggered  PRT  sequences  may  be  achieved  by  use  of 
FIR  (Finite  Impulse  Response  or  "all-zero")  filters  which  operate  on  the  unevenly 
spaced  data.  The  goals  of  this  section  of  the  report  are  to  develop  a  general  theory 
of  these  filters;  to  identify  those  areas  in  which  their  behavior  differs  from  FIR  filters 
that  operate  on  equispaced  data;  and  to  discuss  the  interaction  of  these  filters  with 
weather  parameter  estimation. 

3.1  STAGGERED  PRT  FILTER  STRUCTURE 

We  begin  by  noting  that  the  filtering  of  unevenly  spaced  data  has  received  little 
attention  in  the  engineering  literature.  H.W.  Thomas  and  N.P.  Lutte  treat  this 
topic  analytically  in  two  papers  which  appeared  in  the  early  1970’s 
[Thomas72a.Thomas74a]  .  The  method  of  analysis  employed  in  those  papers  is 
that  of  the  advanced  Z-Transform.  The  approach  taken  here  will  be  somewhat 
different. 

The  staggered  PRT  sampling  scheme  observes  data  values  at  times 

•  ■  ■  ,  —  (c  i  +2c  2 )  ’  —  (c  i  +C  2)^  >  —  ^  2*^r  »  0  >  c  j  t>r  >  (c  j-t-o  2)&r  >  (2c  i+c  2)  ‘  ’ 

where  C]<c2  •  The  sampling  method  is  depicted  in  Figure  3.1  for  the  case 
c  ,=2,0  2=3  (c>r  has  units  of  seconds  throughout  the  following  discussion).  As 
stated  earlier,  with  cx  and  c2  constrained  to  be  relatively  prime  integers,  frequencies 
in  the  range  (  — (  2dr)~1  ,  (  2<Sr)-1]  Hertz  may  be  unambiguously  observed. 

The  output  of  a  causal  FIR  filter  with  N  coefficients  is  a  weighted  sum  of  the 
current  input  value  and  the  N—  1  most  recent  input  values.  In  the  case  of  an 
equispaced  sampling  scheme  with  an  intersample  spacing  of  T  seconds,  the  M th 
most  recent  sample  is  collected  MT  seconds  previous  to  the  current  sample 
( M  =  0,1,  •  •  •  and  M  =  0  corresponds  to  the  current  sample).  In  the  case  of  a 
staggered  sampling  scheme,  however,  the  time  difference  between  the  sampling 
instants  for  the  current  sample  and  the  M  th  previous  sample  is  a  function  of  the 
current  sampling  instant.  Defining 

d0,o=°  >  d0,i=c2  *  d o,2=*'  1+C2  >  do  z=c  i+2 c 2  ,  d 0,4— 2c  1+2C2  ,  d05=2c  i-t-3c2  ,  *  •  * 

d i,o=0  ,  d\fi=c  1  ,  di,2=c1+c2  ,  d1,3=2c1+c2  ,  d1A=2c1+2c2  ,  d1,5=3c1+2c2  ,  •  •  • 

it  may  be  observed  that  for  a  sample  taken  c2&r  seconds  after  the  immediately 
preceding  sample,  the  difference  in  sampling  instants  between  this  sample  and  the 
M  th  preceding  sample  is  do.M^r  seconds.  Similarly,  for  a  sample  taken 
seconds  after  the  immediately  preceding  sample,  the  difference  in  sampling  instants 
between  this  sample  and  the  M  th  preceding  sample  is  dijw  <^r  seconds.  This  is  illus¬ 
trated  for  the  case  cx—2  ,  c2=3  in  Figure  3.2. 
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Figure  3.1:  Staggered  PRT  sampling  scheme 
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The  sainple-to-sample  alternation  of  the  intersample  spacings  between  a  given 
sample  and  its  predecessors  suggests  that  the  use  of  two  sets  of  filter  weights  would 
provide  the  greatest  degree  of  control  in  filtering  time  series  data  collected  in  this 
manner.  One  set  of  weights,  /i0  —{h0l  }/iy\  would  be  selected  to  achieve  a  desired 
response  -  the  meaning  of  response  will  be  defined  subsequently  -  for  samples 
observed  at  times  {d0 ,  ^l/Io1  prior  to  the  current  sample.  This  set  of  weights  would 
be  utilized  when  the  current  sample  is  that  at  time  A  in  Figure  3.2,  and  at  all  other 
times  when  the  current  sample  is  taken  at  p  (Cj-t-c^^  seconds  relative  to  time  A  , 
pel.  The  second  set  of  weights,  hx  —{hXi },n=q1,  would  be  selected  to  achieve  a 
desired  response  for  samples  observed  at  times  Wu^rl/lo1  prior  to  the  current  sam¬ 
ple.  This  set  of  weights  would  be  employed  when  the  current  sample  is  that  at  time 
B  in  Figure  3.2,  and  at  all  other  times  when  the  current  sample  is  taken  at 
p  (cj+Cyl'V  seconds  relative  to  time  B  ,  p  el.  If  x(t)  is  the  continuous  time  function 
being  sampled,  the  filter  is  depicted  schematically  in  Figure  3.3.  where  we  have 
defined  xk  ( A*  <^1  and 


k=-2 

-C2 

k=-l 

0 

k  =0 

c  1 

k  —1 

Ci+C2 

k  =2 

2c  i  +C  2 

k  =3 

We  will  refer  to  the  two  sets  of  filter  weights  h0  and  A lt  employed  as  in  Figure 
3.3  as  a  Staggered  PRT  Filter  of  Length  N ;  and  to  each  set  of  weights  ht  ,  i  =  0,1  , 
as  the  i  th  partial  response  sequence. 
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Figure  3.3:  Staggered  PRT  filter  architecture 


29 


3.2  STAGGERED  PRT  FILTER  RESPONSE 

The  need  for  two  different  partial  response  sequences,  corresponding  to  the  two 
delay  sequences  ,  i  =  0,1  ,  may  be  demonstrated  by  a  simple  example. 

As  a  preface,  we  introduce  the  Z-Transforms  of  the  ht .  The  one-sided  Z-Transform 
of  a  continuous  tune  function  x  (t )  relative  to  some  sampling  period  T  is  defined  as 

CO 

X(2)4Z{x(i)>  =  £  x  (nT  )z  ~n 

n  =0 


If  x{nT)  has  finite  energy,  this  series  converges  everywhere  outside  of  some  circle, 
centered  at  the  origin  in  the  z-plane.  Due  to  the  relative  primality  of  c  j  and  c  2  ,  the 
greatest  common  remainderless  divisor  of  the  two  sampling  periods  c  jAp  and  c2A p  is 
Ap  .  Thus,  the  Z-Transforms  of  the  partial  response  sequences  are  given  by 

C,  (z  )  4 Z{hi }  =  hi M Z~d‘  k  ,  I  =  0,1 

k  =0 


In  the  following  discussions,  the  region  of  convergence  of  the  Ct-  (z )  is  assumed  to 
include  |z  |  =  1. 

By  evaluating  the  Ct  ( z  )  on  the  unit  circle,  relative  to  the  base  sampling  rate  of 
Ap  ,  we  obtain 


Ciiz) 


N-l  'A  t 

S  '*"‘V 

k  =0 


i  -  0,1 


(3.1) 


This  is  the  complex  frequency  response  of  the  partial  response  sequence  /i,  to  a 
sinusoid  at  radian  frequency  u  ,  sampled  at  intervals  of  Ap  .  We  note  that  the  period 
of  Ct  (z  )  ,  as  a  function  of  w  ,  is  the  length  of  the  interval  (—  1  ,  1]  -  the  same 

frequency  interval  which  may  be  unambiguously  observed  for  any  relatively  prime 
integer  values  of  c  r  and  c2  . 

When  h0=hi  (which  corresponds  to  passing  the  staggered  PRT  sequence 
through  a  conventional  tapped  delay  line),  the  relationships  between  C  0  and  C  j  for 
even  values  of  N  differ  from  those  for  odd  values  of  N .  Consider  the  following  case: 
N  —2  ,  c  !=1  ,  c2=  2  and  h0=hl={ 1,1}  .  This  choice  of  parameters  yields 


C0(e'^)  =  2  ,  zeroes  at  wAj-  = 

z 

C  =  l-He/w6r)“1  ,  zeroes  at  uhp  =  ±n 

Thus,  h0  effectively  acts  as  a  notch  filter,  with  notches  at  ±(4Ap)_1  Hz,  while  hx 
functions  as  a  lowpass  filter  with  notches  at  ±(2Ax)_1  Hz.  With  an  input  sequence 
of,  say  { xk  }  =  {sinCwA*  Ap)}  ,  w  =  ~Af 1  ,  it  is  apparent  that  a  modulation  effect,  due 
to  the  mismatch  of  C0  and  C  x  ,  will  be  introduced  into  the  output  sequence  {y*  }  of 
Figure  3.3  . 
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Alternatively,  for  odd  N .  consider  the  case:  A7  =3  ,  c ,  =1  .  c2=2  and 

hn=hi  ={1,1,1}.  We  have 

C0(ey  T)  =  lMeJ  ^r2+(eJ^)-:i 

C,(e7' T)  =  l+(eJJVMcJ^r3  =  (e;"Vc;(eJ'‘,6r) 

In  fact,  for  odd  N  and  a  real  h0=h  lt  it  is  always  the  case  that 

Cx(z)  =z~'MN^){Ci+C2)Cn(z). 

These  two  examples,  the  example  for  N  =2  in  particular,  point  out  the  advan¬ 
tage  of  individual  specifications  for  h0  and  h They  also  show  that  the  effect  of 
alternation  between  h{)  and  /i  1  is  a  consideration  in  selecting  those  coefficients.  As 
the  autocorrelation  properties  of  the  output  sequence  {yk  }  are  of  primary’  interest,  we 
will  subsequently  be  concerned  with  the  effective  power  transfer  fimctions  of  the 
staggered  PRT  filter. 


3.2.1  Average  Power  Transfer  Function 

The  average  power  transfer  function  of  a  staggered  PRT  filter  at  frequency 
u’e(—  -.Sf1  ,  “tvf1]  ,  |G(-A|2  ,  will  be  defined  as  the  filter’s  sample-averaged, 

magnitude-squared,  output  given  an  input  of  a  sampled  phasor  at  radian  frequency 
a1.  Thus  we  have  as  input  {xk  }  =  {eJ '"****},  and 


|GU|2  4  lim 


£  bk  Is 


-  AT-oc  2iVf  +1  k 

For  k  even,  we  note  that  -\._z  =  A*.  —d0,i  and  hence  that 


N- 1 

yt  =  £ 
/=o 


j  \  v,V 


j  \  U’V  ,  -jd-O.l 

e  nO,le 

1=0 


C0{eJ  **) 


Similarly,  for  k  odd,  yk  =  eJ 1(e7w6r)  so  that  equation  3.2  may  be  rewritten  as 

|GM|2=*»  LLwTr|c'(eJ"Sr>|2  +  A  2rar|c^"at,|2(h 

k  odd  k  even 
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|G  (a-)  |2  =  lim 

M  -«-oo 


1  j 

(M  +1)  |C0(e^  ’■  )  |2  +  M  \C  |2 

■  ,  M  even 

2M+1 

1 

M  | C0(ei“h)  |2  -f-  (M+l) |2 

■  ,  M  odd 

2M+1 

Recalling  that  C0(ey  <5r)  and  C1(eJ,"Ar)  are  finite  for  aU  values  of  u:  and  taking  the 
limit,  we  arrive  at  the  expected  result 


|GU)  |2  =  y  |C0(eJuir)  |2  +  j  | Cl(ej"fr)  |2  (3.4a) 

If  the  {hi  j,  }k~i  ,  i  =  0,1  ,  are  real,  we  obtain 


-i-C0fe)Co(-)  +  4-C1fe)C,(i) 

IT ,  f- 

— - ,  - rC 

2  2  2  2 

,  =L.J  '4r 

&p  (*^p 

|G(w)  | 2  represents  the  average  power  gain  of  the  staggered  PRT  filter  for  a 
phasor  input  at  discrete  frequency  jj  e  (—  ~<>f 1  ,  *]  .  Superposition  may  be 
employed  to  determine  the  filter’s  average  power  gain  for  more  complicated  input 
functions. 


3.2.2  Convolutions  of  Staggered  PRT  Sequences 


The  purpose  of  this  section  of  the  report  is  to  examine  the  convolution  operator 
in  the  context  of  Staggered  PRT  sequences.  Following  sections  will  employ  the  con¬ 
volution  operator  in  deriving  their  results. 

Consider  two  equispaced  data  sequences,  {gi ,*  }kJu  >  1  =  0,1  ,  where  the  glk 
are  identically  zero  for  all  k  £{  0,  .  .  .  ,N  —  1 }  .  Denoting  the  convolution  operator  by 
an  inline  *  ,  we  have  by  definition 


fi  =g  o  *  S  i 


N- 1 

2  g  0  ,kg  1 J  —k 
k  =0 


fi  is  then  identically  zero  outside  of  {0,  •  •  ■  ,2(N  —  1)}  and,  subject  to  a  particular 
choice  of  theg(  k  ,  may  be  non-zero  for  any  value  of  l  within  that  set. 

When  two  Staggered  PRT  sequences  are  convolved,  the  time  extent  of  the  con¬ 
volution  is  constrained  by  the  time  extent  of  the  two  sequences,  as  above.  However, 
the  special  structure  of  these  sequences  imposes  additional  constraints  on  the 
instants  at  which  the  convolution  may  be  non-zero.  We  will  demonstrate  this  effect 
for  a  convolution  required  in  the  following  section,  hx  k  *  h 0_k  .  h0  _k  is  the  partial 
response  sequence  h0Je  reflected  about  t=  0.  Recalling  that  5?  is  the  greatest  com¬ 
mon  remainderless  divisor  of  the  two  intersample  delays,  and  c2& r  >  define  the 
following  two  equispaced  sequences  with  intersample  delays  of  c4>. 
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.  i  —0 


Whence, 


,m  Kn  =— c 


h()'k  >  m  —  —d-o.k 
h  ]  j.  ,  ni  =(!])*  ,  i  =1 

0  ,  otherwise 


f  i  =h  1Jk  *  h  o._fc 


OC- 

l  Xj  ^l.Z— m^O,— m 
m  =—oo 


,  /el 


where  /  indexes  sampling  times  in  units  of  &r  .  This  f  z  will  be  identically  zero  for 
values  ofU[-dQtN^,  •  •  •  ,di,w_i]. 

In  particular,  this  ft  may  take  on  non-zero  values  only  at  sampling  times  1 8? 
for  which  le[-d 0jv_i»  •  •  •  ,di,/v-i]  and  one  of  the  foDowing  four  equations  is 
satisfied. 


1)  /  — p  (Cj+C2)  =  0 

2)  /  — p  (C  j  +C  2 1  =  — 0 2 

3)  /  - /)  =  c1  ’  p 

4)  l-p(cx+c 2)  =  C j  — C 2 

Values  of  l  satisfying  these  equations  for  \p  |  <  3  are  displayed  in  Table  3.1  . 


Values  of  1  Satisfying  Eqn.  3.5 

P 

1 

2 

3 

4 

O 

o 

3(c  i+c  2) 

3 c  ^  +2c  2 

4c  J 

+3c  2 

4c  ^  +2c  2 

2 

2(c  j  -fc  2) 

2c,+c2 

3c  | 

+2  c  2 

3c i+c2 

1 

C  i+C2 

Cl 

2c 

1  +c2 

2c  1 

0 

0 

c2 

Cl 

Ci-c2 

-1 

—  (c  j  -tc  2) 

—  C 1  —  2c  2 

c2 

-2c  2 

-2 

—  2(c  j  -fc  2) 

-2^  — 3c2 

-Cl- 

-2c  2 

~  Ci— 3c  2 

-3 

—  3(c  j  -fc  2) 

—  3c  j  —  4c  2 

-  2c  r 

_  3c  2 

—  2cj  — 4c2 

TABLE  3.1 

To  use  Table  3.1  to  determine  those  values  of  l  at  which  f  t  may  be  non-zero, 
we  must  select  a  value  for  N  and  bear  in  mind  that  l  must  lie  in  the  range 
f~c/o,N-l>  ’  ‘  ‘  • 

Two  specific  examples  are  provided  below  for  cx  =  2  ,  c2  =  3  and  N  <4.  For 
this  range  of  N,  p  takes  on  values  of  -1,0,1  . 

Case  1:  N  =  3  ,  f  z  may  be  non-zero  for  (from  Table  3.1) 

(c i+C 2)  C 2  C j  —  C 2  0  C |  2c,  cx+c2 

1  =  -5  -3  -1  0  2  4  5 
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Case  2:  N  =  4  ,  f t  may  be  non-zero  for  (from  Table  3.1) 

-0^202  -2c2  —  ( C  ]  +C  2 )  —  C  2  C,-C2  0  c\  2c  T  C}-hC2  2c  j  +c  2 
1  =  -8  -6  -5  -3  -1  0  2  4  5  7 

Note  that  in  both  Case  1  and  Case  2,  ft  may  be  non-zero  for  values  of  the 
index  /  at  which  h0  _k  and  hlk  are  both  identically  zero,  in  particular  for  l  =  — 1,4 
in  both  cases  and  for  /  =  —  6  in  Case  2.  Also  note  that  Case  2  employed  all  the 
entries  of  Table  3.1  for  p  —  —1,0,1  while  Case  1  did  not. 

As  the  Staggered  PRT  filter  generates  output  values  at  times  A*.  <Sj«  ,  these 
unique  convolution  properties  of  Staggered  PRT  sequences  may  not  be  readily 
apparent  in  the  filter’s  time  domain  output  sequence.  In  the  frequency  domain,  how¬ 
ever,  their  effects  are  evident  -  the  Z-Transforms  of  these  convolutions  will  contain 
powers  ofz  =eJ’'“Sr  other  than  {  ±d(),o  ,  •  •  •  ,  ±d0jv-i  >  ±d i.i  ,  •  '  •  ,  ±di,w_i  )  • 

In  the  following  sections  of  this  report,  we  will  employ  the  notational  conveni¬ 
ence  of  denoting  all  convolutions  by  an  inline  *  .  The  reader  should  bear  in  mind  the 
results  of  this  section  for  convolution  operations  involving  one  or  more  data  sequences 
with  staggered  intersample  delays. 


3.2.3  Lagct  and  La gc2  Power  Transfer  Functions 

Given  some  wide-sense  stationary  input  sequence  {xk}  to  a  linear  time- 
invariant  filter,  the  autocorrelation  function  of  the  output  sequence  {yk  }  is  completely 
characterized  by  the  filter’s  power  transfer  function  and  the  autocorrelation  function 
of  the  input  sequence.  Unfortunately,  |G(cz)  j2  does  not  play  an  exactly  analogous 
role  for  the  Staggered  PRT  filter.  To  understand  why  this  is  so,  we  must  examine 
the  definition  of  the  autocorrelation  function. 


For  some  wide-sense  stationary,  equispaced,  input  sequence,  {xk },  the  auto¬ 
correlation  function  of  {xk  }  at  lag  p  is  defined  to  be  Rxx  (p )  —  Ex  {xkxk+p}  ,  where 
Ex  denotes  an  expectation  with  respect  to  x  and  a  superscript  *  represents  complex 
conjugation.  If  {xk  >  is  passed  through  a  linear  time- invariant  filter  with  impulse 
response  if  i)tLol  and  corresponding  frequency  response  F (uj),  the  autocorrelation 
function  of  the  output  sequence  is  given  by 


Ryy(p)  —Ex 


N  -1 

Ti  xk-lfl 
l  =*0 


1  pv-l 

Zj  xp  +k  -l  f  l 
,/=  0 


=  R, 


fi  *  f-i 


Denoting  the  power  spectra  of  {**}  and  {yfe>  as  Sxx(u.’)  and  Syy  (<U  ,  respectively,  we 
obtain  the  equivalent  result 

SyyU)  =  IF(u’)I2S„(x>) 
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Since  |FU)  |2  is  phaseless,  the  phase  characteristic  of  SVVU’)  is  identical  to  that  of 
S,,U). 

For  the  remainder  of  the  current  section  of  this  report,  we  will  use  k  to  denote 
an  offset  in  sample  indices  which  define  the  lag  at  which  an  autocorrelation  function 
or  cross-correlation  function  is  computed.  Referring  to  Figure  3.3,  we  see  that  for 
even  values  of  k  (  i.e.  for  values  of  A*  which  are  integer  multiples  of  (c1+c2)  ).  the 
autocorrelation  function  of  the  output  sequence  {yk  }  of  a  Staggered  PRT  filter  is 
related  to  the  autocorrelation  function  of  the  input  sequence  {xk }  in  the  following 
manner  (recall  the  preceding  section  of  this  report). 

i  if*'1  H  v-i  FI  [7.v-i  v -i  i* 

Ryyi^f,  )  =  E.V  X  Xn^ll0.t  XI  X*  -M-lho.l  +  X  xn~{h  1/  X  1./ 

(  1=0  n  even  l  =0  i  =0  n  odd 

=  ~  *  hot  /lo.-/)  +  *  hu  *  A l.-/) ]  j*  .  h  even 

or 

Syv  ( =  |  G  ( *■)  1 2SXX  ( J)  ,  for  la gs  p  (c  2  +c  2) ,  P  el 

For  odd  values  of  k  ,  however,  the  power  transfer  functions  which  relate  Svv  (...•) 
and  Sxx(x)  depend  upon  the  particular  value  of  A*.  as  the  autocorrelation  of  iyk  }  is, 
in  these  cases,  actually  a  cross-correlation  of  the  outputs  of  the  two  partial  response 
sequences  h()  and  X-  For  the  case  A*  =p  (c14-c2)+c2  ,  p  el  ,  the  autocorrelation 
function  of  {yk  }  is  given  by 

(  (N- 1  (N-i 

^-vv  1  ~  Ey  1  X  X  %p  (c  ]+c2)+C2+n  -l  ^  0,/ 

l  W=0  (/=0 

=  Rxx  *  hu  *  h q_[  k  ,  k  odd,  A(,  =p(ci+c2)+c2  ,  pel 
or 

Syy(J)  =  c[)(ej  "Sv)C  i{eJ^^)Szx(u)  ,  for  lags  p  (c!+c2)+c2  ,  p  el 

Similarly,  for  the  case  \  =  p  (c  j  +c  2)-t-c :  ,  p  el  ,  the  autocorrelation  function  of  {yt } 
is  given  by 

Ryy  (Aft  )  =RXX  *  hQl  *  h\  _t  k  ,  k  odd  ,  Aft  =p  (c1+c2)+c1  ,  pe I 
or 

Svv  =  C n(ej *fj)C  \  (eJ "'^)SXX  (oJ)  ,  for  lags  p  (c1+c2)+c1  ,  p  el 
Collecting  our  results,  we  have  the  following  for  p  el 


,  k  even  ,  \  =  p(c1+c2) 


Ryv  ( At )  — 


Rr 


Rr 


(R.rx  *  R(),l  ’’’  ^  0.-/)  +  (R.xx  *  h\,l  *  ^  1,-/)  ] 


hu  *  h'0  _i  k  ,  k  odd  ,  A*  =p  (c1  +c2)+c2 


(3.6) 


ho ./  * 


,  k  odd  ,  A*  =p  (c1+C2)+Ci 


Syy  (w) 


J  G  ( w)  1 2StI  (w)  ,  for  lags  p  (c :  +c 2) 
Co(e-/“Sr)C1(e-'  "6r)S„(cj)  ,  for  lags  p  (c1+c2)4c2 
Co(eJ*ST)Cl(eJuST)Sxx(v)  ,  for  lags  p  (Ci+c^+Ci 


(3.7) 


The  fact  that  the  Staggered  PRT  filter  has  the  three  effective  power  transfer 
functions  displayed  in  equation  3.7  is  a  direct  consequence  of  the  fact  that  it  is  not  a 
time-invariant  system.  The  number  of  effective  power  transfer  functions  is  equal  to 
the  number  of  positive  lags  embedded  in  {y*  }  which  are  less  in  magnitude  than  the 
lag  corresponding  to  the  basic  period  of  the  sampling  scheme  -  c1+c2. 

An  analysis  similar  to  that  of  this  section  and  the  preceding  one  could  be  per¬ 
formed  for  a  sampling  scheme  which  possessed  intersample  spacings  of 

>  —  1  *  i  h  »  ^2  ’  *  »  —  1  *  1 1  »  1 

relative  to  some  time  base,  &r  >  where  the  ct  are  relatively  prime.  In  this  case,  the 
analogue  of  the  filter  depicted  in  Figure  3.3  would  be  comprised  of  N  partial 
response  sequences  and  would  possess  N  (N  — 1)+1  effective  power  transfer  func¬ 
tions.  For  the  case  N>2  ,  however,  certain  intersample  delays  could  possess  multi¬ 
ple  effective  power  transfer  functions.  As  an  example,  if  N  =3  and 
Cj=2  ,  c2=3  ,  c3=5  ,  lags  of  c3  and  ^-l-c-  would  both  correspond  to  a  delay  of  56j  ; 
and  the  net  power  transfer  function  for  estimates  of  the  autocorrelation  function  of 
the  filter’s  output  at  a  delay  of  56r  would  be  the  average  of  the  lag  c:+c2  and  lag  c3 
effective  power  transfer  functions  -  if  all  available  sample  pairs  in  {yk}  with  an 
intersample  spacing  of  5^r  were  to  be  used  for  the  estimate. 

For  lags  other  than  p  (c1-fc2),  equation  3.7  shows  that  the  transfer  functions  of 
the  Staggered  PRT  filter  may  have  phase  characteristics  that  are  not  identically 
zero.  In  designing  Staggered  PRT  filters,  then,  it  is  highly  desirable  that  the  magni¬ 
tude  and  phase  characteristics  of  h0  and  hi  be  closely  matched,  so  that  filter’s 
effective  transfer  function  appears  to  be  independent  of  the  particular  lag  at  which 
an  autocorrelation  value  of  the  filter’s  output  sequence  is  to  be  estimated.  The 
interaction  of  the  phase  responses  of  Co(eJ  ^T)C1(eJ  wAr)  and  C0(eju^)C\  (eJ'^)  with 
pulse-pair  velocity  estimation  will  be  discussed  in  section  3.3.2. 
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3.3  PROPERTIES  OF  THE  C,  (2  )  AND  EFFECTIVE  POWER 
TRANSFER  FUNCTIONS 

The  purpose  of  this  section  of  the  report  is  to  provide  a  qualitative  discussion  of 
some  of  the  frequency-domain  properties  of  the  Staggered  PRT  filters. 


3.3.1  Z-Plane  Zero  Distribution  of  |G(u.)  |2 

When  the  hl  k  are  real,  the  z-plane  zeroes  of  |GM|2  (more  precisely,  the 
zeroes  of  C 0(z  )C 0(z_1)  +C  x{z  )C  1(2~1)  )  display  both  conjugate  symmetry  and 
reciprocal  symmetry  about  the  unit  circle.  To  see  this,  note  that  real  h,  j.  yield 
Ct-  (2 )  with  zeroes  occurring  in  conjugate  pairs.  The  zeroes  of  C,  (z-1)  will  be  the 
inverses  of  those  of  C,-(z),  hence  the  zeroes  of  Ciiz  JC^z-1)  -  the  totality  of  the 
zeroes  of  Ct(z)  and  C;  (2  -  will  possess  conjugate  and  reciprocal  symmetry. 

Zeroes  of  |G(~)  |2  will  correspond  to  one  of  two  cases.  Either 
C0(z  )C0(z~1)  =C  1(z  )C  1(z~1)  =  0  or  C 0(z  )C0(z~l)  =  —C  x(z  )C  x(z  -1).  In  both  cases, 
conjugate  and  reciprocal  symmetry  follow  from  the  properties  of  the 
Ci(z)Ci{z~l)  ,i  =0.1.  These  symmetries  result  in  a  |G(w)|2  which  is  a  real, 
even,  fimction  of  u.  Examples  of  these  z-plane  zero  distributions  will  be  presented 
subsequently. 


3.3.2  Phase  Responses  of  (e7"^)C  1(e7utSr)  and  C0(eJ^^)C \  (e7*^) 

As  mentioned  earlier,  the  filter’s  effective  power  transfer  functions  for  lags  other 
than  p  (c  l+c2)  ,  p  el  may  have  non-zero  phase  characteristics.  The  goal  of  the  next 
several  paragraphs  is  to  describe  conditions  on  the  selections  of  stagger  ratio  and 
filter  length  which  are  sufficient  to  guarantee  linear  phase  characteristics  for 
Co  (e7w5r)C  ife7^)  and  C^ie^' T)C  *  (eJ  '^)  .  In  the  present  context,  we  need  only 
discuss  one  of  these  transfer  functions,  C  q  (e7  ^)C  ](e7  ''"’r)  ,  as  phase  linearity  of  one 
guarantees  phase  linearity  of  the  other.  A  linear  phase  characteristic  is  desirable, 
as  its  only  effect  is  to  induce  a  known  time  shift  into  the  autocorrelation  sequence  at 
the  filter’s  output. 

Referring  to  equations  3.6  and  3.7,  we  see  that 


Cm0(ejuh)C1(ejuh)  =Z{hu  *  h’o^} 


j. 


This  is  the  lag  c2  effective  power  transfer  fimction,  and  it  will  possess  a  linear  phase 
characteristic  if  fm  —hii  *  ho-i  is  a  time  function  which  is  real  and  even  with 
respect  to  some  time  instant. 

If  we  examine  the  sampling  instants,  m  6j- ,  mel  ,  at  which  f  m  may  be  non¬ 
zero  (  e.g.  via  Table  3.1);  and  then  successively  difference  those  sampling  times  to 
obtain  the  intersample  delays  in  fm ,  we  see  the  following  pattern  (differences  are 
displayed  left  to  right  from  the  most  negative  pair  of  sampling  instants  at  which  fm 
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may  be  non-zero  to  the  most  positive  pair). 
N=2: 


C j  ,  C2— C  i  ,  C x 

N=3: 

C1  ,  Cv  ,  C2—C  1  ,  Cj  ,  Cj  ,  C2-C] 

N=4: 

C  i  ,  C2—C  i  ,  CX  ,  Ci  ,  C2—C  i  ,  C  j  ,  C2—C  i  ,  C  i 

N=5: 

c  1  ,  Cj  ,  C  2— c  1  ,  C]  ,  Cl  ,  C2—C  !  ,  C  i  ,  C  i  ,  C  2—C  \  ,C]  ,  Cj  ,  C  2~C  \ 

N=6: 

C  i  ,  c  2— C  i  ,  C  i  ,  C  i  ,  C2- C  i  ,  C  i  ,  C  i  ,  c2— ■ C  j  ,  C  i  ,  C  i  ,  C2—C  ]  ,  C  i  ,  C  ]  ,  c  2 — C  i  C  i 

For  values  of  N  <6  ,  fm  only  possesses  even  intersample  delays  for  even  values 
of  N  .  We  conjecture  that  fm  possesses  even  intersample  delays  for  all  even  values 
of  N  and  non-even  intersample  delays  for  all  odd  values  of  N  .  For  N  =  2,4,6  the 
central  intersample  delay  of  c2— Ci  is  centered  at  time  —  (c2— c1)6r/2  ,  which  implies 
the  possibility  of  a  linear  phase  characteristic  of  exp  (+j^(c2— ^ c1)6r/2)  for 
Cq  (eJ^)C  .  A  linear  phase  characteristic  will  be  obtained  if  the  values  of 

fm  possess  even  symmetry  about  the  central  intersample  delay. 

We  note  that  the  convolution  of  two  sequences  with  odd  symmetry  always 
results  in  a  sequence  with  even  symmetry.  This  may  be  easily  demonstrated.  Con¬ 
sider  two,  finite  energy,  equispaced  data  sequences  {gti*  >  1  =  0,1  which  pos¬ 

sess  odd  symmetry  -  gik  =  —gi  _k  ,  i  =0,1.  Their  convolution,  evaluated  at 
some  index  l  is  given  by 

oo  oo 

So,k  *  8i,k  i  =  I]  8\,k8o,i-k  =  I]  8\,-k8oj;~i 

k  =— oo  k  =— oo 

—oo 

=  X)  8i,k8o,-i-k  =8o,k  *  g\,k  _/ 

k  =+co 

As  will  be  seen,  it  is  often  the  case  that  the  design  of  high- pass  Staggered  PRT 
filters  with  narrow  stopbands  (  e.g.  clutter  filters)  results  in  partial  response 
sequences  which  possess  odd  symmetry  for  even  values  of  N  ,  in  the  sense  that 

hi,o  =  —hijj- 1  ,  hiA  =  —hi^r_ 2  ,  •  •  •  ,  i  =  0,1 

Combining  the  results  of  the  preceding  several  paragraphs,  it  is  apparent  for  even 
values  of  N  and  partial  response  sequences  with  odd  symmetry  that  the  lag  c  }  and 
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lag  c2  effective  transfer  functions  -  C0(eJ"  '1)C  ^  (eJ~  ')  and  C (eJ'  r)C  ^eJ  ‘  T)  -  will 
possess  linear  phase  characteristics  of  exp  (—j  ( c 2  — c  1 ) at /  2 )  and 
exp  (+Ju.  (c2— Ci)t>r/  2)  ,  respectively. 

Figures  3.4  through  3.7  display  these  phase  properties  for  N  =3  and  A’  =4  with 
c  !=1  and  c2=  2  .  The  particular  values  for  Ci  and  c2  were  selected  for  the  purpose  of 
demonstration,  only.  Figures  3.4  and  3.5  display  |G(cj)  |2  and  Cn(eJ'“*r)C  i  (eJ^)  , 
respectively,  for  the  case  N  =3  ;  Figures  3.6  and  3.7  provide  the  corresponding  plots 
for  the  case  N  = 4  (and  partial  response  sequences  with  odd  symmetry).  The  format 
for  these  plots  is  as  follows.  The  abscissa  is  the  positive  portion  of  the  Nyquist  inter¬ 
val  in  Doppler  velocity  space,  for  a  3  GHz  (S-Band)  radar  with  =  0.45455  ms. 
The  corresponding  Nyquist  velocity  is  55  m/s.  The  left-hand  ordinate  axis  is  power 
gain  in  dB,  while  the  right-hand  ordinate  axis  in  Figures  3.5  and  3.7  is  phase  shift 
in  degrees.  The  block  of  text  below  the  plot  on  the  left  displays  the  stagger  ratio, 
radar  constants,  <Sr  and  N .  The  block  of  text  below  the  plot  on  the  right  displays  the 
parameters  used  in  designing  the  partial  response  sequences  -  design  techniques  are 
the  topic  of  subsequent  sections  of  this  report.  Figures  3.4  and  3.6  also  are  marked 
to  indicate  the  Nyquist  velocities  corresponding  to  PRT’s  of  c1^r,c2t'kr  and  . 
Those  markings  are  Vua>1  ,  Vua/,  and  VuaAlrlfo!ded  ,  respectively.  As  c  t=l  for  these 
examples,  the  labels  Vua>1  and  Vuaunfo|ded  are  overlaid. 

Comparing  Figures  3.4  and  3.5,  we  see  that  the  gain  characteristics  of 
|GU’)  |2  and  C()(eJU'r)C  *  (eJ,"<'T)  are  well  matched.  This  is  desirable  as  it  implies 
that  the  Staggered  PRT  filter’s  gain  effects  on  the  power  spectrum  of  its  output  will 
be  essentially  independent  of  the  particular  output  value,  or  collection  of  values, 
examined.  As  expected,  however,  the  phase  response  of  C  0{e^^)C\(eJ  is  non¬ 
linear.  The  non-linearity  is  especially  benign  for  this  example  -  more  realistic  values 
of  c1:c2  (  e.g.  2:3,  3:4,  3:5,  5:7,  etc.)  tend  to  exacerbate  the  non-linearity. 

In  comparing  Figures  3.6  and  3.7,  the  gain  characteristics  of  \G(ui)\2  and 
C0(eJ"Sr)C  j  (eJ~  T)  may  be  observed  to  differ.  For  even  values  of  N  ,  this  disparity 
can  become  acute  for  more  realistic  values  of  c}:c2  .  The  phase  response  of 
C0(eJufr)C  *  (e7^)  is  linear,  with  the  expected  slope  of  —  90(c2-c1)  degrees  per  k 
radians  of  digital  frequency. 

Figures  3.8  through  3.10  provide  an  example  of  velocity  estimation  from  the 
autocorrelation  sequence  of  the  staggered  PRT  filter’s  output.  De-aliasing  is  accom¬ 
plished  as  follows: 

1)  a  lookup  table  containing  all  combinations  of  rate  and  c2&r 
velocity  folds  for  velocities  in  the  range  ±uu  =  ±A(4<5r)_1  is  precom¬ 
puted  "offline"  -  Table  3.2  is  an  example  of  the  lookup  table  for 
c !  =2  ,  c2=3  and  v  U=1  m/s  (see  Figure  1.4,  also); 
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'52963  ’2 


(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 
FILTER  ID:  123 


STAGGER  RATIO  1 .2  AUTOCORRELATION-BASED  DESIGN: 

8t  IS  0.45455  ms  INPUT  S/C  POWER  RATIO:  -1 0.0  dB 

PRF  1  =  2200.0  Hz,  PRF  2  =  1 100  0  Hz  OUTPUT  S/C  POWER  RATIOS:  71 .6  dB. 

3  FILTER  COEFFICIENTS  71.6dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH.  0  00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  ,  00000 

Figure  3.4:  |G(w)  |2  for  a  staggered  PRT  filter,  N=3 
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MAGNITUDE  (dB) 


■5Z963  '3 


MAGNITUDE  AND  PHASE  FUNCTIONS  FOR  C0C  \  vs  DOPPLER  VELOCITY  (m,s) 

FILTER  ID  123 


STAGGER  RATIO  1  2 
6t  IS  0  45455  ms 

PRF  1  =  2200.0  Hz.  PRF  2  =  1 1 00.0  Hz 

3  FILTER  COEFFICIENTS 

RADAR  FREQUENCY  IS  3.00000  GHz 


AUTOCORRELATION-BASED  DESIGN 
INPUT  S/C  POWER  RATIO:  -10.0  o3 
OUTPUT  S  C  POWER  RATIOS:  71  6  dB. 

71  6  dB 

NORMALIZED  CLUTTER  WIDTH  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS 
0.00000  0.00200  C.  85000  1.00000 


Figure  3  5:  Magnitude  and  phase  functions  for  C  0C  * ,  N  =3 
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PHASE  (deg) 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  124 

STAGGER  RATIO  1 :2  AUTOCORRELATION-BASED  DESIGN: 

ST  IS  0.45455  ms  INPUT  S'C  POWER  RATIO:  -10.0  dB 

PRF  1  =  2200  0  Hz.  PRF  2  =  1100  0  Hz  OUTPUT  SC  POWER  RATIOS:  109  3  dB. 

■  FILTER  COEFFICIENTS  1 1 0.5  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0  00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 

Figure  3.6:  | G  ( co)  | 2  for  a  staggered  PRT  filter.  N  =4 
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MAGNITUDE  (dB) 
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MAI  ,NI  llltH  AND  PHASI  I  I  INC;  I  IONS  H)RC„C  ",  v-,  DOPP1  F  H  Vf  IOCIIY  (m  ■'.) 

f  ii  r f  m  a)  cm 

SI  ACiCit  M  FFATIO  1  :*  AUTOOORRFl  A1  ION  BASF  L)  DF  SIGN 

1-1,  IS  0  4 04 S')  ms  INPUT  S  O  POWFH  RATIO  10  0  dB 

I’RI  1  -  FS’OO  (I  M/.  f’RI  -  1 100  0  H/  ()l  IFPI  IT  S  C  POWI  R  HA  NOS  100  I  dB. 

4  I  II  TFR  COI  I  I  1011  NTS  1 10  S  dB 

RADAR  F  HFOUFNGY  IS  1  00000  GH/  NOFIMAl  l/F D  01  UTTFR  WIDTH  0  00.'V:i 

NORMA!  I/F.DW  SF’I  OTRIJM  BRE  AKPOIN I ! 

o  ooooo  o  oo;>oo  o  hsooo  i  ooooo 


PHASE  (deg) 


Folding  Table:  C]=2  ,  c?=3  and  v„=l  m/s 


Rate  c  j  Xp 

Number  of  Folds 

1 

0 

0 

0 

1 

Correction  (m/s) 

-1 

0 

0 

0 

+1 

Rate  c2Xp 

Number  of  Folds 

1 

1 

0 

1 

1 

Correction  (m/s) 

-sy.3 

-2/3 

0 

2/3 

2/3 

TABLE  3.2 

2)  potentially  aliased  lag  Cj  and  lag  c2  velocity  estimates,  Uj  and  v2 
respectively,  are  generated  from  the  filter’s  output  sequence; 

3)  using  the  table  computed  in  (1),  the  magnitude  of  the  difference  in 
each  combination  of  corrected  estimates,  |  ulc  —  v2c  I,  is  computed; 
and 

4)  the  final  velocity  estimate  is  the  average  of  the  ulc  and  v2c  which 
minimized  |u  lc  —  v2c  j  in  (3). 

Figure  3.8a  displays  |G(^')  j"  for  a  Staggered  PRT  filter  with  ^=1,  c2=2  and 
N  =2.  v  u  is  55  m/s.  The  filter  was  designed  using  the  OSCRM  design  technique 
which  will  be  described  in  section  4.3.  Figure  3.8b  displays  the  lag  cx  transfer  func¬ 
tion  C  n(eJ^)C\  (eJ  “  T).  N  is  even  and  the  /i,  are  odd,  in  the  sense  that  hit0=-hit lt 
hence  C0(eJutiT)C\  (eJ"T)  possesses  a  linear  phase  response.  Figure  3.9a  displays 
the  uncorrected  lag  c  t  (solid  line)  and  lag  e2  (dashed  line)  velocity  estimates  versus 
true  velocity  -  the  weather  spectrum  width  is  1  m/s  for  this  example.  These  were 
computed  via  a  equation  1.4c  as  — X(4ne16rrWg[i?;v;v  (c^)]  and 
— X(4/7C2Xr)_1arg[/i’vv  (c26r)],  respectively.  Properly  aliased  velocities  would  have 
resulted  in  line  segments  with  slopes  of  unity  in  this  figure.  Instead,  the  slopes  are 
V 2  and  fi/2. 

Noting  that  the  phase  responses  of  C0(eJ^)C  *  (e'/u<'r)  and  Cq  (ey“^)C  i(e^&T) 
correspond  to  an  advance  and  delay  of  14^,  respectively,  we  may  instead  compute 
the  lag  Cj  and  lag  c2  velocity  estimates  as 

~  4»~-^ar8l^lcAI1  ,3'8a) 

=  4^WSr‘rgl^’,C‘ir>1  (3'8a) 

to  account  for  the  offset  in  Ryy  relative  to  Rxx .  The  resulting  estimates  are  plotted 
in  Figirre  3.9b.  These  estimates  are  properly  aliased,  with  the  exception  that  the  lag 
c2  estimate  folds  at  approximately  22  m/s,  rather  than  at  V4uu  =  27.5  m/s.  The 
cause  of  the  premature  fold  is  currently  a  mystery;  as  it  occurs  in  a  velocity  region 
where  |  Cn(eJ^)C  \  (eJU&T)  |  is  not  constant,  a  partial  explanation  may  lie  in 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  122 

STAGGER  RATIO  1 :2  AUTOCORRELATION-BASED  DESIGN: 

Sr  IS  0.45455  ms  INPUT  S  C  POWER  RATIO:  -10.0  dB 

PRF  1  =  2200  0  Hz.  PRF  2  =  1 1 00.0  Hz  OUTPUT  S/C  POWER  RATIOS:  28.6  dB. 

2  FILTER  COEFFICIENTS  34.0  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1  00000 


Figure  3.8a:  |G(vV)  \2  for  an  OSCRM  design,  ,  c2=  2  and  N  —2 
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MAGNITUDE  AND  PHASE  FUNCTIONS  FOR  C0C;  vs  DOPPLER  VELOCITY  (m,s) 

FILTER  ID:  122 

STAGGER  RATIO  1 :2  AUTOCORRELATION-BASED  DESIGN: 

5t  IS  0.45455  ms  INPUT  S/C  POWER  RATIO:  -1 0  0  dB 

PRF  1  =  2200.0  Hz,  PRF  2  =  1 100.0  Hz  OUTPUT  S/C  POWER  RATIOS:  28.6  dB, 

2  FILTER  COEFFICIENTS  34  0  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 


Figure  3.8b:  Lag  ca  transfer  function  for  OSCRM  design  of  figure  3.8a 
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PHASE  (deg) 


'52963  -6 


VELOCITY  ESTIMATES  (m/s)  vs  DOPPLER  VELOCITY  (m/s) 
FILTER  ID:  122 


WEATHER  SPECTRUM  WIDTH  IS  1 .0  m/s 
LAG  1  RESET  TO  1.0 
LAG  2  RESET  TO  2.0 
STAGGER  RATIO  1:2 
ST  IS  0.45455  ms 

PRF  1  =  2200.0  Hz.  PRF  2  =  1 100.0  Hz 

2  FILTER  COEFFICIENTS 

RADAR  FREQUENCY  IS  3.00000  GHz 


Figure  3.9a:  Uncorrected  lag  Cj  and  c2 

<rl  =1  m/s 


NYQUIST  SPECTRAL  FOLDING  ENFORCED 
AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10.0  dB 
OUTPUT  S/C  POWER  RATIOS:  28  6  dB. 

34.0  dB 

NORMALIZED  CLUTTER  WIDTH:  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 


velocity  estimates  from  filter  of  figure  3.8, 
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VELOCITY  ESTIMATES  (m/s)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  122 

WEATHER  SPECTRUM  WIDTH  IS  1 .0  m/s  NYQUIST  SPECTRAL  FOLDING  ENFORCED 
LAG  1  RESET  TO  0.5  AUTOCORRELATION-BASED  DESIGN: 

LAG  2  RESET  TO  2.5  INPUT  S/C  POWER  RATIO:  -10  0  dB 

STAGGER  RATIO  1 :2  OUTPUT  S/C  POWER  RATIOS:  28.6  dB. 

8T  IS  0.45455  ms  34.0  dB 

PRF  1  =  2200.0  Hz.  PRF  2  *  1 100.0  Hz  NORMALIZED  CLUTTER  WIDTH.  0.00273 

2  FILTER  COEFFICIENTS  NORMALIZED  W-SPECTRUM  BREAKPOINTS: 

RADAR  FREQUENCY  IS  3.00000  GHz  0.00000  0.00200  0.85000  1 .00000 


Figure  3.9b:  Corrected  lag  c x  and  c2  velocity  estimates  from  filter  of  figure  3.8, 

av  =1  m/s 
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MEAN  VELOCITY  ESTIMATE  ERROR  (m/s)  vs  MEAN  WEATHER  VELOCITY  (m/s) 

FILTER  ID:  122 


WEATHER  SPECTRUM  WIDTH  IS  1.0  m/s. 

LAG  1  RESET  TO  0  5 

LAG  2  RESET  TO  2.5 

STAGGER  RATIO  1  2 

?>t  IS  0  45455  ms 

PRF  1  =  2200.0  Hz.  PRF  2  =  1 100.0  Hz 
MINIMUM  DISPLAYED 
VELOCITY  2  0  m/s 


2  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 
NYQUIST  SPECTRAL  FOLDING  ENFORCED 
AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10  0  dB 
OUTPUT  S/C  POWER  RATIOS:  28.6  dB. 

34.0  dB 

NORMALIZED  CLUTTER  WIDTH:  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS 
0.00000  0.00200  0.85000  1.00000 


VELOCITY  DE-ALIASING  OPERATING  UPON  LAG  1  VELOCITY  ESTIMATES  AT  EACH  PULSE  RATE. 

Figure  3.10:  Errors  in  final  de-aliased  estimates  generated  from  lag  c,  and  c2 
estimates  of  figuie  3.bb,  ov  -1  m/s 
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distortion  of  the  weather  spectrum.  The  effect  of  the  premature  fold  is  a  positive 
bias  of  approximately  11  m/s  (2(27.5-22))  for  the  aliased  lag  c2  estimate.  Velocity 
de-aliasing  performs  properly,  however  -  the  errors  in  the  final  de-aliased  estimates 
are  displayed  in  Figure  3.10.  Note  that  the  averaging  performed  in  step  (4)  reduces 
the  bias  in  the  unfolded  estimate  to  5.5  m/s  for  velocities  greater  than  22  m/s.  A 
more  intelligent  de-aliasing  algorithm,  with  knowledge  of  the  lag  c  x  and  lag  c2 
transfer  functions,  could  further  reduce  that  bias. 

For  those  C0(eJI‘,'t)Ci(e‘,'''r)  which  meet  all  of  the  aforementioned  criteria  for 
linear  phase,  but  possesses  zeroes  with  odd  multiplicity  on  the  unit  circle,  180° 
discontinuities  will  be  introduced  into  the  otherwise  linear  phase  characteristic.  In 
such  cases,  equations  3.8  will  not  yield  properly  aliased  estimates.  Figure  3.11  pro¬ 
vides  an  example  of  a  C0(e;“^)C]  (e;‘‘<'r)  with  a  discontinuous,  "linear",  phase 
characteristic.  This  is  an  OSCRM  design  with  2,  c2=3  and  N  =4.  The  unambi¬ 
guous  velocity  is  55  m/s  and  the  weather  spectrum  width  is  1  m/s.  Figure  3.12 
displays  the  lag  C\  and  lag  c2  velocity  estimates  of  equations  3.8.  The  effects  of  the 
phase  discontinuity  in  C 0(eJ'~''iT)C  \  (eJ'~^)  at  «34  m/s  is  clearly  evident  -  note  that 
correctly  aliased  estimates  are  obtained  for  velocities  below  the  phase  discontinuity, 
although  the  lag  c2  estimate  folds  prematurely  at  »15  m/s  rather  than  at  I8V3 
m/s.  The  appropriate  folding  table  for  this  stagger  ratio  and  unambiguous  velocity 
is  Table  3.2  with  the  entries  labeled  Correction  (m/s)  all  scaled  by  55.  Applying  the 
de-aliasing  algorithm  to  these  estimates,  we  obtain  the  final  estimate  errors  shown 
in  Figure  3.13.  The  errors  in  the  lag  cx  and  lag  c2  estimates  of  Figure  3.12  result 
in  the  selection  of  an  improper  additive  correction  (-55  m/s  in  this  case)  for  velocities 
in  the  range  of  34  to  55  m/s.  Empirically,  discontinuous  "linear"  phase  characteris¬ 
tics  like  that  of  Figure  3.11  are  more  common  than  truly  linear  phase  characteristics 
like  that  of  Figure  3.8b. 

This  section  of  the  report  may  be  summarized  as  follows.  For  Staggered  PRT 
clutter  filters,  even  values  of  N  generally  result  in  linear  phase  characteristics  for 
the  lag  Cj  and  lag  c2  effective  power  transfer  functions  and  disparities  in  the  gain 
characteristics  of  |  G  ( F)  |2  and  C 0(eJ"^)C  *  (e7^)  .  When  the  phase  characteristics 
are  truly  linear,  equations  3.8  will  yield  properly  aliased  lag  c  x  and  lag  c2  velocity 
estimates.  Odd  values  of  N  generally  result  in  non-linear  phase  characteristics  for 
the  lag  cl  and  lag  c2  effective  power  transfer  functions  and  good  agreement  among 
the  gain  characteristics  of  |G  (w)  |2  and  C 0(eJ"&r)C *  (eJu^)  -  as  mentioned  in  section 
3.2,  odd  N  and  a  real  hn=hx  guarantee  that  f  C0(ejw&T)  |=|  C\(e*uSr)  j  and  hence 
that  |G  (w)  |2=  |  C u(eJJ>T)C  1  (e1  jSr)  |  .  All  other  things  held  constant,  increasing  the 
values  of  N  ,Cj  or  c2  tend  to  exacerbate  phase  non-linearities  and  magnitude 
response  mismatches  in  a  filter  design. 
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MAGNITUDE  AND  PHASE  FUNCTIONS  FOR  C0C',  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  234 

STAGGER  RATIO  2:3  AUTOCORRELATION-BASED  DESIGN: 

6t  IS  0.45455  ms  INPUT  S/C  POWER  RATIO:  -10.0  dB 

PRF  1  =  1 100 .0  Hz.  PRF  2  =  733.3  Hz  OUTPUT  S/C  POWER  RATIOS:  95.8  dB. 

4  FILTER  COEFFICIENTS  96.1  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH.  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0  85000  1.00000 


Figure  3.11:  Lag  c  j  transfer  function  for  OSCRM  design,  cT=2  ,  c2=3  and  N  =4 
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PHASE  (deg) 
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VELOCITY  ESTIMATES  (m/s)  vs  DOPPLER  VELOCITY  (m/s) 
FILTER  ID:  234 


WEATHER  SPECTRUM  WIDTH  IS  1.0  m/s 

LAG  2  RESET  TO  1.5 

LAG  3  RESET  TO  3.5 

STAGGER  RATIO  2:3 

fcT  IS  0.45455  ms 

PRF  1  =  1 100.0  Hz.  PRF  2  =  733.3  Hz 
4  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 


NYQUIST  SPECTRAL  FOLDING  ENFORCED 
AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10  0  dB 
OUTPUT  S/C  POWER  RATIOS:  95  8  dB. 

96.1  dB 

NORMALIZED  CLUTTER  WIDTH:  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS. 
0  00000  0.00200  0.85000  1.00000 


Figure  3.12:  "Corrected"  lag  C!  and  c2  velocity  estimates  from  filter  of  figure  3.11, 

aL  =1  m/s 
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MEAN  VELOCITY  ESTIMATE  ERROR  (m/s)  vs  MEAN  WEATHER  VELOCITY  (m/s) 

FILTER  ID:  234 

WEATHER  SPECTRUM  WIDTH  IS  1 .0  m/s  4  FILTER  COEFFICIENTS 
VELOCITY  DE-ALIASING  OPERATING  RADAR  FREQUENCY  IS  3.00000  GHz 
UPON  LAG  1  VELOCITY  ESTIMATES  NYQUIST  SPECTRAL  FOLDING  ENFORCED 
AT  EACH  PULSE  RATE  AUTOCORRELATION-BASED  DESIGN 

LAG  2  RESET  TO  1 .5  INPUT  S/C  POWER  RATIO:  -10.0  dB 

LAG  3  RESET  TO  3.5  OUTPUT  S/C  POWER  RATIOS:  95  8  dB. 

MINIMUM  DISPLAYED  VELOCITY:  2.0  m/s  96.1  dB 

STAGGER  RATIO  2:3  NORMALIZED  CLUTTER  WIDTH.  0.00273 

5t  IS  0  45455  ms  NORMALIZED  W-SPECTRUM  BREAKPOINTS. 

PRF  1  =  1100.0  Hz.  PRF  2  =  733.3  Hz  0.00000  0.00200  0.85000  1.00000 


Figure  3.13:  Errors  in  final  de-aliased  estimates  generated  from  lag  Cj  and  c2 
estimates  of  figure  3.12,  av  =1  m/s 
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3.3.3  Uncontrollable  Zeroes  in  the  Ct  (z ) 

In  contrast  to  the  design  of  FIR  filters  operating  on  equispaced  data,  not  all  of 
the  Staggered  PRT  filter’s  z-plane  zeroes  may  be  directly  specified.  The  number  of 
z-plane  zeroes  for  an  equispaced  FIR  filter  is  equal  to  the  number  of  coefficients  it 
possesses,  excluding  the  coefficient  at  lag  zero  which  is  a  scalar  multiplicative  factor, 
as  shown  below. 

A(z )  =a0(  1  +a1z-1  + a +  ■  •  •  +  aN_lzN~1) 

In  selecting  the  coefficients  for  such  a  filter,  then,  the  designer  may  directly  specify 
the  location  of  each  z-plane  zero.  A  staggered  PRT  filter,  however,  generally 
possesses  many  more  zeroes  than  filter  coefficients.  With  a  stagger  of  3:5  (  i.e. 
c i=3, c 2=5  ),  for  example,  the  C,  (z  )  have  the  form  displayed  below. 

C 0(z  )  =hi)  {)  +h0  lz  5  +  hQ  2z  8  d"  ^ 0,3-  13  +  ‘  ‘  '  +^o.a,-i2  °"V  ' 

C  1(2  )  =  h\  o  +  h  1 '  \z  ^  •+■  h  ^  2 z  8  +  h  j  3 z  ' '  +  •••  +  1  _  1  z  1 1 

For  every  two  coefficients  added  to  the  partial  response  sequences  beyond  h0  0  or  h  i  0 , 
an  additional  Cj+c2  z-plane  zeroes  are  added  to  C0(z  )  or  C  j(z  )  ,  respectively,  only 
two  of  which  are  controllable.  We  note  that  N  must  be  odd  if  C0(z  )  and  C  1(z )  are 
to  ha,_e  the  same  number  of  zeroes.  If  N  is  odd,  the  C,  (z )  are  both  of  order 

VAN  —1)  (c  j  -Hz  2 )  -  If  N  is  even,  however,  C0(z)  is  of  order  c2+((N  /2)~  l)(c  i+c2) 

while  Cj(z  )  is  of  order  Cj-KflV/2)- l)(c  i+c2)  • 

The  stagger  ratio,  itself,  is  an  important  parameter  of  a  Staggered  PRT  filter 
design  for  two  reasons. 

1)  The  stagger  ratio  determines  the  number  of  uncontrollable  zeroes 
incurred  by  partial  response  sequences  for  a  given  value  of  N  . 

2)  For  filter  designs  in  which  the  controllable  zeroes  are  located  in  the 
z-plane  near  z  =1  (  e.g.  for  clutter  filters  or  other  high-pass  filters 
with  narrow  stopbands),  the  stagger  ratio  appears  to  determine  the 
z-plane  distribution  of  the  uncontrollable  zeroes.  These  are  the 
zeroes  which  determine  the  filter’s  passband  characteristics,  such  as 
ripple. 

Figures  3.14  through  3.16  provide  examples  of  the  effect  of  stagger  ratio  choice 
on  the  z-plane  zero  distribution  of  |G(w)  |2  .  The  partial  response  sequences  were  all 
designed  to  place  the  controllable  zeroes  of  the  C,  (z  )  at  z  =1.  Each  figure  consists  of 
two  plots,  "a"  and  "b".  The  ”a"  plots  are  z-plane  zero  plots  of  |G  (uO  |2  -  the  unit  cir¬ 
cle  is  displayed  on  these  plots,  along  with  a  block  of  text  on  the  lower  left  describing 
the  basic  filter  parameters.  The  "b"  plots  are  corresponding  plots  of  the  magnitude 
response  of  |GM|2  ,  with  the  same  plot  format  as  Figure  3.4  .  Each  "b"  plot 
displays  |G(u/)  |2  for  an  entire  Nyquist  interval  -  digital  frequencies  from  0  to  2 n  - 
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although  different  values  of  dj.  were  employed  in  generating  the  plots.  N  =3  for  all 
three  figures,  implying  two  controllable  zeros  and  (c i+c2)—  2  imcontrollable  zeros  for 
each  partial  response  sequence.  For  these  examples,  the  ht  k  are  real. 

Figure  3.14  displays  | G  ( -  )  |2  for  the  case  Cj=2  ,  c2=3  .  With  N  odd,  C„(z  ) 
and  C  i(z  )  possess  the  same  numbers  of  intersample  delays  of  Cj  and  c2  and  hence 
the  same  distributions  of  powers  of  z  in  the  products  C0(z)C0(z_1)  and 
C^zlCjlz-1)  .  Thus,  |G(a)|2  also  possesses  this  distribution  of  powers  of  z  , 
which  extends  from  z  iri+c,i)  to  zCl+C2  .  For  the  particular  values  of  c,  and  c2  ,  the 
distribution  extends  from  z~5  to  z5  ,  therefore  |G(cc)  |2  possesses  10  zeroes. 

Referring  to  Figure  3.14a,  we  see  that  4  of  these  zeroes  are  located  at  z—1. 
The  controllable  zeroes  of  C0(z )  and  C \{z)  ,  which  were  all  placed  at  z  =1  ,  carry 
over  to  | G  ( -h  |2  for  the  following  reason.  The  C,fz  )  are  of  the  form 

4.V-1-2 

C1(z)=(z-1)2  ri  .  i  =0,1 

m  =  1 


Thus 


1 

C,(z  )Ct  ( — )  =  (z  — 1  )2(z  _I— l)2  n  (z -zt  m  )(z_1-zI>m  )  .  i  =0,1 

2  m  =1 

and 

I G  ( -  )  | 2  =  ~C  (/z  )C0(— )  +  ICjlzlCili) 

Z  z  Z  z 


2 

(z  — l)2(z_1— l)2  £ 

n  iz -zlym )(z-i-z,,« ) 

i=i 

m  =1  J 

The  remaining  six  zeroes  in  Figure  3.14a  are  the  imcontrollable  zeroes  resulting 
from  the  z,  m  in  the  preceding  equation  and  give  rise  to  the  dips  in  the  curve  of  Fig¬ 
ure  3.14b  at  35.  52  and  70  m / s  Doppler  velocity. 

Figure  3.15  displays  |G(c.>  |2  for  the  case  cx=3  ,  c2=4.  As  in  Figure  3.14, 
N  =  3  and  all  controllable  zeroes  of  the  partial  response  sequences  are  placed  at 
z  =1  .  In  this  case,  however,  |G(co)  |2  possesses  10  imcontrollable  zeroes.  This  fact 
is  reflected  in  the  increased  number  of  passband  ripples  in  Figure  3.15b  relative  to 
Figure  3.14b. 

Figure  3.16  displays  |G(w)  |2  for  the  case  ^=5  ,  c2=6.  All  other  parameters 
of  the  filter  design  are  identical  to  those  of  the  preceding  two  figures.  |G(u>)|2 
possesses  18  imcontrollable  zeroes  in  this  case. 

Increasing  the  values  of  ct  and  c2  does  not  necessarily  worsen  the  passband 
characteristics  of  a  Staggered  PRT  filter.  The  effect  of  the  imcontrollable  zeroes  is  a 
function  of  both  „heir  number  and  their  z-plane  locations.  Certain  stagger  ratios  ( 
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Z  PLANE  ZERO  PLOT  OF  q>  (z)G,  d'z)+C,  (z)C,  (1  /z)  FOR  FILTER  233 
NYQUIST  FREQUENCY  IS  (26T  )'’Hz 

STAGGER  RATIO  2:3  4  ZEROES  AT  z  =  1 

3  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 


Figure  3.14a:  Z-Plane  zeroes  of  |GU>)  |2  -  c,=2  ,  c2=3  ,  AT  =3  and  all  controllable 
zeroes  at  z  -1 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  35 

STAGGER  RATIO  2:3  3  FILTER  COEFFICIENTS 

ftT  IS  047619  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  1050.0  Hz.  PRF  2  =  700.0  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  3.14b:  | G(*v)  |2  -c:='2  ,  c2= 3  ,  N  =3  and  all  controllable  zeroes  atz  =1 
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2-PLANE  ZERO  PLOT  OF  C0(z)C0(1/z)+C,  <z)C,  (1/Z)  FOR  FILTER  63 
NYQUIST  FREQUENCY  IS  (28  T)  _1Hz 

STAGGER  RATIO  3:4  4  ZEROES  AT  z=1 

3  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.15a:  Z-Plane  zeroes  of  |G(^’)  |2  -  Cj=3  ,  c2= 4  ,  N  =  3  and  all  controllable 
zeroes  at  z  —1 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  37 

STAGGER  RATIO  3:4  3  FILTER  COEFFICIENTS 

&T  IS  0.35714  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  933.3  Hz.  PRF  2  =  700.0  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  3.15b:  \G  (  J)  \2  -  c}=3  ,  c2=4  ,  N  =3  and  all  controllable  zeroes  at  z  =1 
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Z-PLANE  ZERO  PLOT  OF  C0(z)C0(1/Z)+C,  (z)C,  (1/z)  FOR  FILTER  68 
NVQUIST  FREQUENCY  IS  (2ST  ) Hz 

STAGGER  RATIO  5:6  4  ZEROES  AT  z=1 

3  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.16a:  Z-Plane  zeroes  of  |G(w)  |2  -  Cj=5  ,  c2=&  ,  AT=3  and  all  controllable 
zeroes  at  z  =1 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 
FILTER  ID:  38 


STAGGER  RATIO  5:6 
8T  IS  0  23810  ms 

PRF  1  =  840.0  Hz.  PRF  2  =  700.0  Hz 


3  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 
MAXIMALLY  FLAT  DESIGN 


Figure  3.16b:  |G(u>) 


c1=5  ,  0  2=6  ,  N= 3  and  all  controllable  zeroes  at  z  =1 
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e.g.  5:7  )  result  m  passband  characteristic  with  small  ripples  and  a  steep  slope  in 
the  transition  region  between  the  stopband  and  passband.  Placing  all  controllable 
zeroes  at  z  1  makes  the  transition  region  slope  a  function  of  the  uncontrollable 
zeroes,  as  well. 

Empirically,  the  selection  of  a  staler  ratio  appears  to  determine  the  /.-plane 
loci  for  |(/  (u.)  |2’s  uncontrollable  zeroes.  Within  the  unit  circle  in  the  upper  half- 
plane,  the  loci  appear  to  lie  second  order  functions  such  as  hyperbolas  and  parabolas. 
Outside  of  the  unit  circle  and  in  the  lower  half-plane,  the  loci  are  determined  by  the 
conjugate  and  reciprocal  symmetries  discussed  earlier.  To  date,  an  analytical 
description  of  these  loci  as  a  function  of  the  stagger  ratio  has  not  been  obtained. 
Factoring  polynomials  is  generally  difficult. 

As  an  example  of  such  loci,  Figures  3.17a  through  3.17f  display  the  /.-plane 
positions  assumed  by  uncontrollable  zeroes  which  are  successively  added  to  \(H  |  . 

Figures  3.l8a-f  are  corresponding  plots  of  | ( »(...>)  | 1  versus  Doppler  velocity.  Each 
design  displayed  in  these  figures  possesses  a  5:7  stagger  ratio  with  all  controllable 
zeroes  clustered  at  z  1  -  this  type  of  design  will  he  discussed  in  section  4.2  as  the 
Maximally  Flat  liesjxmttc  design  technique.  The  values  of  N  for  Figures  .'5.17a  fare 
3,  4,  5,  (),  7  and  H,  respectively;  resulting  in  4,  (5,  8,  10,  12  and  14  controllable 
zeroes,  and  20,  .‘52,  40,  52,  (50  and  72  uncontrollable  zeroes,  respectively.  The  radi- 
als  in  the  upper  half-plane  of  the  zero  plots  are  included  for  the  purpose  of  illustra¬ 
tion,  only.  As  N  increases,  additional  uncontrollable  zeroes  appear  on  the  loci, 
increasing  the  depth  of  the  passband  ripples  of  |(»  (u)  p.  The  ripple  locations  are 
essentially  fixed.  The  additional  zeroes  placed  near  z  1  increase  the  .stopband 
width  with  N ,  and  the  transition  region  width  (e.g.  the  velocity  interval  over  which 
|G(u4|2  goes  from  -(»()  <114  to  -3  <114 )  decreases,  at  a  lesser  rate,  with  N.  As  the 
stopband  width  increases  more  rapidly  than  the  transition  hand  width  decreases,  the 
passband  edge  also  increases  with  increasing  N . 
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Z-PLANE  ZERO  PLOT  OF  C0(z)C0(1/z)+C,  (z)C,  (1/z)  FOR  FILTER  573 
NYQUIST  FREQUENCY  IS  (25T  )  Hz 

STAGGER  RATIO  5:7  4  ZEROES  AT  z  =  1 

3  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.17a:  Z-Plane  zeroes  of  |GM  |2  -  c-[=5  ,  c2= 7  ,  N  =3  and  all  controllable 
zeroes  at  z  =1 
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2-PLANE  ZERO  PLOT  OF  C0(z)C0(1/z)+C,  (z)C,  (1/2)  FOR  FILTER  574 
NVQUIST  FREQUENCY  IS  (26T  ) *’  Hz 

STAGGER  RATIO  5:7  6  ZEROES  AT z  =  1 

4  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.17b:  Z-Plane  zeroes  of  |G(cj)  |2  -  c}=5  ,  c2=7  ,  iV=4  and  all  controllable 
zeroes  at  z  =1 
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2  PLANE  ZERO  PLOT  OF  C0(Z)C0(1/z>+C,  (z)C,  (1/2)  FOR  FILTER  575 
NYQUIST  FREQUENCY  IS  (2  5T) Hz 

STAGGER  RATIO  5  7  2  ZEROES  AT  z  =  1 

5  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.17c:  Z-Plane  zeroes  of  |G(u?)  |2  -  c^—5  .  c2=7  ,  N  —  5  and  all  controllable 
zeroes  at  z  =1 
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Z-PLANE  ZERO  PLOT  OF  C0  (z)C0(1/Z)+Ct  (z)C,  (1/2)  FOR  FILTER  576 
NYQUIST  FREQUENCY  IS  (26T  ) Hz 

STAGGER  RATIO  5:7  0  ZEROES  AT  z  =  1 

6  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.17d:  Z-Plane  zeroes  of  |G(u)  j2  -  Cj=5  ,  c<i=l  ,  N  —6  and  all  controllable 
zeroes  at  z  =1 
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Z-PLANE  ZERO  PLOT  OF  CQ  (z)C 0(1/z)^C1  (2)0,  (1/2)  FOR  FILTER  577 
NYQUIST  FREQUENCY  IS  (25T  )  *’Hz 

STAGGER  RATIO  5:7  0  ZEROES  AT  z  *  1 

7  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3.17e:  Z-Plane  zeroes  of  |G(u,)  |2  -  C}=5  ,  c2= 7  .  N  =7  and  all  controllable 
zeroes  at  z  =1 
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Z-PLANE  ZERO  PLOT  OF  CQ(z)CQ{Vz)+C,  (z)C,  (1/z)  FOR  FILTER  578 
NYQUIST  FREQUENCY  IS  (2ST  ) Hz 

STAGGER  RATIO  5:7  0  ZEROES  AT  z  =  1 

8  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 

Figure  3 . 1 7 f:  Z-Plane  zeroes  of  |G(cu)  \2  -  c,=  5  ,  c2=  7  ,  N=  8  and  all  controllable 
zeroes  at  z  =1 
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(MAGNITUDE)  (dB)  vs  DOPPLER  VELOCITY  (m  s) 

FILTER  ID:  573 

STAGGER  RATIO  5:7  3  FILTER  COEFFICIENTS 

St  IS  0.15455  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  440.0  Hz.  PRF  2  =  314.3  Hz  MAXIMALLY  FLAT  DESIGN 


|G(^.)  !"  versus  Doppler  velocity  -  c,=5  ,  c2=:7  ,  N  =3  and  all  con 
trollable  zeroes  at  z  =1 
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STAGGER  RATIO  5:7  4  FIL  ter  COEFFICIENTS 

f>r  IS  0  45455  ms  RADAR  FREQUENCY  IS  3  00000  GHz 

PRF  1  =  440  0  Hz.  PRF  2  =  314  3  Hz  MAXIMALLY  FLAT  DESIGN 


'itfure  3.18b: 


|G(w)  | 2  versus  Doppler  velocity  -  c 
trollable  zeroes  at  z  =1 


5  ,  c-,-7  ,  N  =4  and  all  con 
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(MAGNITUDE)2  <dB)  vs  DOPPLER  VELOCITY  (m  s) 

FILTER  ID:  575 

STAGGER  RATIO  5:7  5  FILTER  COEFFICIENTS 

ftT  IS  0  45455  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  440. L  Hz.  PRF  2  =  314  3  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  3.18c:  |G(w)  |2  versus  Doppler  velocity  -  Cj=5  ,  c-i—l  .  N  =5  and  all  con¬ 

trollable  zeroes  at  z  =1 
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10  20  30  40  50  60 


(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m  s) 

FILTER  ID:  576 

STAGGER  RATIO  5:7  6  FILTER  COEFFICIENTS 

§T  IS  0.45455  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  440.0  Hz.  PRF  2  =  314.3  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  3.18d:  | GU>  |2  versus  Doppler  velocity  -  c,=5  ,  c2=7  ,  N  =6  and  all  con¬ 

trollable  zeroes  at  z  =1 
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(MAGNITUDE)  (dB)  vs  DOPPLER  VELOCITY  (m  s) 

FILTER  ID:  577 

STAGGER  RATIO  5:7  7  FILTER  COEFFICIENTS 

f>T  IS  0  45455  ms  RADAR  FREQUENCY  IS  3  00000  GHz 

PRF  1  =  440.0  Hz.  PRF  2  =  314.3  Hz  MAXIMALLY  FLAT  DESIGN 

Figure  3.18e:  |G(lu)  |2  versus  Doppler  velocity  -  Ct=5  ,  c2=7  ,  Ar=7  and  all  con 
trollable  zeroes  at  z  =1 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID.  578 

STAGGER  RATIO  5:7  8  FILTER  COEFFICIENTS 

&T  IS  0  45455  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  440  0  Hz.  PRF  2  =  314.3  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  3.18f:  |G(~')  |2  versus  Doppler  velocity  -  Cj=  5  ,  c2=  7  ,  N  =8  and  all  con¬ 

trollable  zeroes  at  z  =1 
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4.0  STAGGERED  PRT  FILTER  DESIGN  TECHNIQUES 

The  current  section  of  this  report  will  describe  several  techniques  for  determin¬ 
ing  the  staggered  PRT  filter’s  partial  response  sequences.  Although,  the  techniques 
will  be  discussed  in  the  context  of  clutter  filter  design,  the  Coefficient  Selection  and 
OSCRM  techniques  have  general  applicability. 

4.1  DESIGN  TECHNIQUE:  COEFFICIENT  SELECTION 

A  wide  variety  of  tools  exist  for  the  design  of  FIR  filters  which  operate  on 
equispaced  data,  most  notably  the  Parks-McClellan  algorithm  which  was  developed 
in  the  early  1970’s  [Parks72a]  .  The  Coefficient  Selection  technique  allows  these 
equispaced  FIR  design  tools  to  be  employed  in  the  design  of  Staggered  PRT  filters. 
Coefficient  Selection  operates  as  follows. 

1)  Values  are  selected  for  ,  c  j  and  c2  ,  and  for  the  Staggered  PRT 
filter  length,  N  . 

2)  Equispaced  FIR  filter(s)  are  designed  with  the  desired  response  for 
an  intersample  delay  (PRT)  of  Ap  .  If  N  is  odd,  a  single  FIR  design 
with  di  i+l  coefficients  (  i.e.  an  impulse  response  duration  of 
di,N-i^r  seconds)  is  required.  If  N  is  even,  matched  FIR  designs 
with  do,Af-i+l  and  d  i  ,v_j  +1  coefficients  (impulse  response  dura¬ 
tions  of  d0jv_i&r  and  dltN_-ldT  seconds,  respectively)  are  required. 

3)  The  Staggered  PRT  filter’s  partial  response  sequences  are  deter¬ 
mined  by  time-aligning  those  sequences  with  the  FIR  design(s)  and 
selecting  coefficients  from  the  FIR  design(s)  at  the  matching  sam¬ 
pling  instants.  Denoting  the  FIR  design  for  the  case  where  N  is 

odd  as  {fihLt',  we  have 

*o,o  =  f  o  -  *o,i  =  f doA  ,  •  •  •  ,  ho.N-i  =  f dos_ , 

*  i,o  =  fo  >  *1,1  =  fdt !  *  *'■  >  *  l.W-1  =  fdiy_x 

In  the  case  that  N  is  even,  and  denoting  the  FIR  designs  with 
impulse  response  durations  of  do,y-i*r  and  duy-i^r  as  if  i°h  =o  -1 
and  {f  z1}/ Jgv-1  ,  respectively,  we  have 

*0.0  ~f  O'  *o,i  =  f  d0i  >  '  '  '  >  *o.w-i  =  f  doy_  j 

*1,0  =  fo  >*1,1  =  /'d11  >  •••  1  *1JV-1  = 

Given  Cj  ,  c2  ,  N  and  Ap,  the  staggered  PRT  partial  response  sequence  gen¬ 
erated  from  an  equispaced  FIR  design  by  Coefficient  Selection  possesses  the  Discrete 
Fourier  Transform  (DFT)  which  most  closely  approximates  that  of  the  equispaced 


FIR  design  in  the  sense  of  the  l2  norm.  More  concisely,  the  DFT  is  an  isometric 
mapping  from  /2[0 ,M  —1]  onto  itself. 

l2[0,M  —  1]  is  the  Hilbert  Space  of  comp] ex -valued  sequences  .  For 

any  two  elements  of  /2[0,M-1],  rl  and  r2  ,  the  /2[0,M  —  1]  inner  product  and  norm 
are  defined  as 

M- 1 

(r1  |r2)  4  •£  n'rf* 
l=o 

jjr1  ||  —  4/(r!  |r]) 


respectively.  It  is  clear  that  the  DFT  maps  /2[0,M  —  1]  to  itself. 

Define  the  DFT  and  Inverse  DFT  (IDFT),  respectively,  for  /2[0,M  —1]  as 

(  ,  M-l 

DFT {rt }  ±{Rk}  =  V  rie-j2M/M 

I  VM  i  _q 


IDFT  {Rk  }  4  }  = 


1 

VM 


M- 1  I 

V  Rkei 

*  =o  J 


To  establish  that  the  DFT  is  an  isometric  mapping  of  Z2[0,M  —  1]  to  itself,  we  must 
show  that  |jr1-r2||  =  ||  for  all  r \r2  el2[0M — 1]  -  Now, 
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j 27, k  {I  -m  )/M 

l_(ej2*k(l-m  )/M  ^ 


=  0  ,  TUtU 


Therefore, 

HR1-#2  |j2  =  Ijr1— r2  ||2  (4.1) 

Since  r1  and  r2  were  arbitrary  elements  of  Z2[0,M-1] 
DFT:Z2[0,M-1]— *72[0,M  —  1]  is  proven  to  be  an  isometry.  This  is  a  slightly  modified 
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statement  of  ParsevaFs  Theorem. 

For  a  given  stagger  ratio  and  Staggered  PRT  filter  length,  equation  4.1 
guarantees  that  each  of  the  partial  response  sequences  generated  by  the  Coefficient 
Selection  technique  most  closely  approximate  the  transfer  fimction(s)  of  the 
corresponding  equispaced  FIR  design(s)  at  the  frequency -domain  sampling  points  of 
the  DFT,  in  the  sense  of  the  /2[0,M—  1]  norm.  Equation  4.1  does  not,  however,  pro¬ 
vide  an  indication  of  the  Staggered  PRT  filter’s  characteristics  at  digital  frequencies 
other  than  {  2rk  / M  ,  nor  does  it  provide  guidance  in  the  choice  of  stagger 

ratio. 

The  desired  frequency  characteristic  of  clutter  filters  for  weather  radar  applica¬ 
tions  is  that  of  a  high-pass  filter  with  a  narrow  stopband,  relative  to  the  Nyquist 
interval,  and  a  short  transition  region  between  the  stopband  and  passband. 
Equispaced  FIR  filters,  with  some  odd  number  M  of  coefficients,  designed  for  these 
applications  typically  possess  impulse  responses  such  as  the  one  depicted  in  Figure 
4.1.  in  which 


ft  = 


M  —  1 
o  ,  l 


l  =0,  •  •  •  ,M— 1 

M  —1 
2 


0  ,  otherwise 


l  9^ 


M  — 1 
2 


where  o  is  a  multiplicative  constant.  When  an  equispaced  FIR  filter  of  this  sort  is 
employed  for  a  Coefficient  Selection  design  of  a  Staggered  PRT  filter,  some  guidance 
in  the  selection  of  a  stagger  ratio  is  provided  by  the  fact  that  the  zeroes  of  a  complex 
polynomial  are  a  continuous  function  of  the  polynomial’s  coefficients  [Marden49a]  . 

Consider  two  equispaced  impulse  response  sequences  of  length  M  ,  and 

{ gi  }/l,7 1  ,  and  define  their  Z-transfonns  to  be  F(z  )  and  Glz  )  ,  respectively.  Let 

ti  =gi  -ft  ,  l  =  0,  •  •  •  ,M  -1 


Then  the  magnitude  of  the  difference  in  their  Z-transforms  may  be  bounded  as  fol¬ 
lows 


\F  (z)-G(z) 


M- 1 

£ 


£iZ' 


M -1  Af-1 

<  £  M"'  I  <  fmax  £  I*”'  I 


f=0 


/=0 


/=0 


where  cmax  =max{  |>  .  If  we  evaluate  the  Z-transforms  on  the  unit  circle,  z  =eju, 

we  obtain  the  following  upper  bound. 

|F  (ic)— G  (cj)  |  <M £max  (4.2) 

Once  a  value  for  N  has  been  selected,  equation  4.2  provides  the  following  guidance 
for  choosing  a  stagger  ratio.  The  stagger  ratio  should  be  selected  so  that  the 
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Figure  4.1:  Equispaced  FIR  clutter  filter  impulse  response,  M  odd 
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coefficient  of  maximum  modulus  in  the  equispaced  design  is  incorporated  in  both  of 
the  partial  response  sequences.  In  particular,  for  the  equispaced  designs  described  in 
the  preceding  paragraph,  the  central  coefficient  in  the  impulse  response  should  be 
incorporated  in  both  of  the  partial  response  sequences. 

Figures  4.2  through  4.4  provide  a  graphical  display  of  equation  4.2  in  the  case 
of  a  Coefficient  Selection  design.  Each  figure  consists  of  two  plots  "a"  and  "b”.  The 
"a"  plots  are  z-plane  zero  plots  of  |G(ta)  |2  for  the  various  designs;  the  ”b"  plots  are 
the  corresponding  frequency  domain  plots  as  a  function  of  digital  frequency.  As  the 
|G(u;)  | 2  arc  even  functions  of  frequency,  each  "b"  plot  displays  only  the  positive  por¬ 
tion  of  the  Nyquist  interval  <ae[0,7r]. 

Figure  4.2  displays  the  equispaced  FIR  design  employed,  which  possesses  17 
filter  coefficients.  The  filter  was  designed  with  the  Parks-McClellan  algorithm  -  the 
parameters  provided  to  that  algorithm  appear  on  the  bottom  left-hand  side  of  the 
plots.  Each  zero  plotted  in  Figure  4.2a  is  actually  a  pair  of  zeroes.  The  effect  of  the 
zero  pairs  near  r  =1  is  clearly  evident  in  Figure  4.2b;  they  generate  the  sharp  notch 
at  0.03  radians  per  sample  interval.  The  impulse  response  of  this  filter  is  similar  to 
that  shown  in  Figure  4.1. 

Figure  4.3  displays  |G(*r)  |2  for  the  corresponding  Staggered  PRT  design  with 
c  j=3  ,  c 2=5  and  N  -5.  With  this  choice  of  Staggered  PRT  parameters,  both  partial 
response  sequences  incorporate  the  central  FIR  coefficient  (as  h0  2  and  hx2)-  Com¬ 
paring  Figures  4.2a  and  4.3a,  we  see  zero  pairs  for  the  FIR  design  beginning  to 
diverge  in  the  Staggered  PRT  design.  The  divergence  is  relatively  benign  for  all  the 
zero  pairs  except  for  the  two  pairs  near  z  =1  .  The  corresponding  effect  in  Figure 
4.3b  is  that  the  passband  characteristic  of  the  Staggered  PRT  filter  is  close  to  that 
of  the  equispaced  filter,  while  the  depth  of  the  stopband  has  been  considerably 
reduced. 

Figure  4.4  displays  |  G  ( v)  |2  for  the  corresponding  Staggered  PRT  design  with 
Cj=2  .  c2=4  and  N=  6  (note  that  cx  and  c2  are  not  relatively  prime  -  this  is 
reflected  in  the  fact  that  Figure  4.4b  displays  a  period  of  tt  for  |G(u.’)  |2  ,  indicating 
that  this  design  is  equivalent  to  a  1:2  stagger  with  a  base  sampling  rate  of  2<Sr  ). 
In  this  case,  neither  of  the  partial  response  sequences  incorporates  the  central  FIR 
coefficient.  The  divergence  of  the  zeroes  in  Figure  4.4a  relative  to  the  equispaced 
design  is  so  severe  that  the  zero  distribution  has  almost  become  angularly  uniform. 
Figure  <*.4b  depicts  a  filter  which  possesses  an  allpass  characteristic. 

Unfortunately,  the  Z-Transforms  of  partial  response  sequences  generated  by 
the  Coefficient  Selection  technique  possess  a  local  maximum  for  their  error  relative  to 
the  corresponding  equispaced  design  in  the  frequency  domain,  at  ^>=0  -  the  center  of 
the  desired  stopband  for  clutter  filters.  This  was  evident  in  Figure  4.3. 

To  see  why  this  is  so,  let  fi  and  gi  be  as  before,  and  denote  their  Z- 
Transforms,  evaluated  on  the  unit  circle  z  =e>"  ,  by  F  (J)  and  G(  J)  .  We  impose  the 
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ZPLANE  ZERO  PLOT  OF  C0(z)C0(1/Z)+C,  (z)C,  (1/2)  FOR  FILTER  1117 
NYQUIST  FREQUENCY  IS  (28  T) Hz 

STAGGER  RATIO  1:1  M.M  S.E  DESIGN  FROM  FIR 

1 7  FILTER  COEFFICIENTS  NFILT  -  1 7  NBANDS  -  2 

4  ZEROES  ATz=1  EDGE  0.00000  0  00909  0.02727  0  50000 

8t=  1000  00000  ms  WEIGHTS:  10  00  1  00 

Figure  4. ‘2a:  Z-Plane  zeroes  of  squared-magnitude  response  for  17  point 

equispaced  FIR  filter 
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0  0.5  1.0  1.5  2.0  2.5  3.0 


(MAGNITUDE) 2  (dB)  vs  DIGITAL  FREQUENCY  (Radians  Per  Sample) 

FILTER  ID:  1117 

STAGGER  RATIO  1 :1  M.M.S.E  DESIGN  FROM  FIR. 

8t  IS  10000.00000  ms  NFILT  =  1 7  NBANDS  =  2 

PRF  1  =  10Hz.  PRF2=  10Hz  EDGES: 

17  FILTER  COEFFICIENTS  0  00000  0  00909  0.02727  0  50000 

WEIGHTS: 

10  0  1.00 


Figure  4.2b:  Squared-Magnitude  response  for  17  point  equispaced  FIR  filter 
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Z-PLANE  ZERO  PLOT  OF  C0(z)C0(1/z)+C,  (z) C,  (1/z)  FOR  FILTER  355 
NYQUIST  FREQUENCY  IS  <26T  ) Hz 


STAGGER  RATIO  3  5 

M  M  S  E  DESIGN  FROM  FIR 

5  FILTER  COEFFICIENTS 

NFILT  -  1 7  NBANDS  -  2 

0  ZEROES  AT  z  =  1 

EDGE:  0.00000  0  00909  0  02727 

0.50000 

5t  =  1000  00000  ms 

WEIGHTS:  10  00  1.00 

Figure  4.3a:  Z-Plane  zeroes 

of  |  G  (  j)  | "  for  coefficient 

selection  design 

c  j  =3  ,  c  2=0  and  N  =5 
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(MAGNITUDE)2  (dB)  vs  DIGITAL  FREQUENCY  (Radians  Per  Sample) 
FILTER  ID:  355 

STAGGER  RATIO  3:5  M  M  S  E  DESIGN  FROM  FIR: 

6t  IS  10000  00000  ms  NFILT  =  17  NBANDS  =  2 

PRF  1  =03  Hz.  PRF  2  =  02  Hz  EDGES: 

5  FILTER  COEFFICIENTS  0  00000  0  00909  0  02727  0  50000 

WEIGHTS: 

10  0  1  00 


Figure  4.3b: 


|  G  ( w)  | 2  for  coefficient  selection  design  -  c ,  =3 


Co=5  and  N  =5 
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Z-PLANE  ZERO  PLOT  OF  C0(2)C0(1/2)+C,(2)C,  (1/z)  FOR  FILTER  246 
NYQUIST  FREQUENCY  IS  (2  8  r ) Hz 

STAGGER  RATIO  2:4  M.M.S.E  OESIGN  FROM  FIR 

6  FILTER  COEFFICIENTS  NFILT  ■  17  NBANDS  -  2 

0  ZEROES  ATz=1  EDGE:  0.00000  0  00909  0  02727  0.50000 

ST=  1000  00000  ms  WEIGHTS:  10  00  1.00 

Figure  4.4a:  Z-Plane  zeroes  of  |G(cc)|2  for  coefficient  selection  design 
c  ,  —2  .0-2=4  and  N  =6 
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(MAGNITUDE)2  (dB)  vs  DIGITAL  FREQUENCY  (Radians  Per  Sample) 
FILTER  ID:  246 

STAGGER  RATIO  2  4 
ST  IS  •'0000  00000  ms 
PRF  1  =  0  5  Hz.  PRF  2  =  0.3  hz 
6  FILTER  COEFFICIENTS 


Figure  4.4b:  \G  (u>)  |2  for  coefficient  selection  design  -  c  x—2  ,  c2= 4  and  N  =6 


M  M  S.E  DESIGN  FROM  FIR 
NFILT  =  17  NBANDS  =  2 
EDGES: 

0.00000  0  00909  0  02727  0  50000 
WEIGHTS: 

10.0  1.00 
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additional  requirement  that  the  ft  and  gt  are  real  for  all  values  of  l.  Define 
D  {  J)  =  |F(a;)— ( G(*j)  |2  .  D(u:)  is  then  a  real-valued,  even,  function  of  us  and  neces¬ 
sarily  possesses  a  local  extremum  at  w  =  0  .  When  the  equispaced  filter  is  of  the 
type  depicted  in  Figure  4.1,  D  (J)  can  be  differentiated  twice  to  show  that  this 
extremum  is  a  local  maximum. 

Thus,  we  have  the  result  that  an  application  of  the  Coefficient  Selection  tech¬ 
nique,  in  which  the  equispaced  FIR  coefficients  not  selected  are  all  of  the  same  sign, 
will  lead  to  local  maxima  in  the  errors  of  the  Z-transforms  of  the  partial  response 
sequences,  relative  to  the  equispaced  design,  at  cj=0  .  This  limits  the  applicability  of 
the  Coefficient  Selection  technique  in  designing  Staggered  PRT  clutter  filters. 
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4.2  DESIGN  TECHNIQUE:  MAXIMALLY  FLAT  RESPONSE 

The  Maximally  Flat  Response  (MFR)  design  technique  assigns  ail  controllable 
zeroes  of  the  Staggered  PRT  partial  response  sequences  to  z  =1  .  resulting  in  a 
high-pass  design  with  a  deep  notch  at  ^=0  .  This  technique  is  discussed  in  some 
detail  by  Thomas  and  Lutte  [Thomas74a]  . 

The  average  power  transfer  function  |G(a>)|2  of  a  Staggered  PRT  filter  of 
length  N  possesses  2(N  —  1)  controllable  zeroes.  If  all  these  zeroes  are  to  be  placed 
at  z  =1,  then  all  the  controllable  zeroes  of  the  partial  response  sequences  must  also 
be  placed  at  z—1.  The  MFR  technique  achieves  this  distribution  for  the  partial 
response  sequences  in  a  manner  analogous  to  a  Butterworth  filter  design. 

MFR  coefficients  are  independently  generated  for  the  two  partial  response 
sequences  by  setting  C,(z  )  and  its  first  N  —  2  derivatives  to  zero  at  z  =1  for  i  =0,1; 
and  then  solving  the  resulting  system  of  linear  equations  for  ](  •  •  •  . 

Since  the  effect  of  ht  0  in  the  Ct  (z  )  is  that  of  a  multiplicative  constant  (  i.e.  it  does 
not  effect  the  zero  locations  of  C,  (z )  ),  it  may  be  chosen  arbitrarily.  We  will  select 
hi' o=l  ,  i  =0,1.  This  results  in  the  two  systems  of  equations 


where 


1  1 

$.1,1  $',1,2 


$,N-2.1  $,N-  2,2 


1 

$,l.Af-l 


I.N-2JV-1 


hi,\ 

-1 

hi, 2 

0 

hi  jv  -l 

•  o 

,  i  =0,1 


(4.3a) 


'A ,/  ,m 


(4.3b) 


Note  that  the  factorial  growth  in  | $ii/n  |  with  the  row  index,  /,  may  cause  the 
coefficient  matrix  in  4.3a  to  become  ill-conditioned.  This  ill-conditioning  is  evident  in 
the  designs  of  Figures  3.11  where  the  controllable  zeroes  are  clustered  about,  rather 
than  collocated  at,  z  =1  for  N> 3. 

As  MFR  places  all  controllable  zeroes  of  |G(o/)  |2  for  a  Staggered  PRT  filter  at 
z  =1  ,  the  passband  characteristics  of  such  filters  are  wholly  determined  by  the 
stagger  ratio  and  filter  length.  Figure  4.5  displays  |G  (u)  |2  for  an  MFR  design  with 
Cj=4  ,  c2~5  and  N= 3  .  Figure  4.5a  displays  the  z-plane  zero  plot  for  this  filter. 
The  cluster  of  controllable  zeroes  may  be  observed  at  z  =1  .  Figure  4.5b  displays 
|G(-j)  | 2  in  the  Doppler  velocity  domain,  for  a  radar  frequency  of  3  GHz  and  base 
sampling  interval  of  0.26471  ms  (Nyquist  velocity  of  94  m/s),  over  an  entire 
Nyquist  interval.  The  passband  characteristic  of  this  filter  is  a  typical  one  for 
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Figure  4.5a: 


Z-PLANE  ZERO  PLOT  OF  C0(z)q,  (1/z)+C,  (z)C,  (1/Z)  FOR  FILTER  66 
NVQUIST  FREQUENCY  IS  (2ST  ) Hz 

STAGGER  RATIO  4:5  4  ZEROES  AT  Z  =  1 

3  FILTER  COEFFICIENTS  MAXIMALLY  FLAT  DESIGN 


Z-Plane  zeroes  of  \G(u)  |2  for  MFR  design  -  c  j=4  ,  c2=5  and  N  =6 
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0  20  40  60  80  100  120  140  160  180  200 


(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  27 

STAGGER  RATIO  4:5  3  FILTER  COEFFICIENTS 

5t  is  0.26471  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  944  4  Hz.  PRF  2  =  755.5  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  4.5b:  | G(J)  j2  for  MFR  design  -  c  j=4  .  c2=5  and  N  =6 
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01  234567  89  10  11 


(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID.  27 

STAGGER  RATIO  4:5  3  FILTER  COEFFICIENTS 

8t  IS  0.26471  ms  RADAR  FREQUENCY  IS  3.00000  GHz 

PRF  1  =  944  4  Hz.  PRF  2  =  755.5  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  4.5c:  for  MFR  design  ,  c2— 5  and  N  =£ 


Staggered  PRT  filters  with  small  values  of  N .  Figure  4.5c  is  a  magnified  version  of 
Figure  4.5b,  for  Nyquist  velocities  up  to  10  m/s.  Since  scan-modulated,  distributed 
clutter  for  S-Band  weather  radars  has  a  spectrum  width  on  the  order  of  0.25  m/s, 
we  see  that  this  filter  provides  adequate  clutter  rejection.  The  filter  has  a  wide  tran¬ 
sition  region,  however,  with  the  3  dB  passband  edge  located  at  approximately  15 
m/s  or  at  0.16  of  the  Nyquist  velocity.  Even  if  &j>  were  reduced  so  that  the  Nyquist 
velocity  were,  say,  50  m/s,  the  filter’s  stopband  and  transition  band  would  still 
cause  unacceptable  attenuation  of  weather  spectra. 

Figure  4.6  displays  a  Doppler  velocity  domain  plot  for  an  MFR  design  which 
possesses  the  same  stagger  ratio  as  the  design  of  Figure  4.5,  but  with  iV=10.  Note 
the  increased  width  of  the  stopband,  and  the  increased  passband  ripple  resulting 
from  the  additional  uncontrollable  zeroes.  Recall  the  examples  at  the  end  of  section 
3.3.3  in  which  stopband  width  increases  were  achieved  at  the  expense  of  moving  the 
passband  edge  upwards.  For  larger  values  of  N ,  the  first  passband  ripple  -  whose 
location  remains  essentially  unchanged  -  serves  to  fix  the  passband  edge,  but  the 
passband  ripples  are  unacceptably  deep  in  these  cases. 

The  MFR  design  technique  is  of  limited  utility  in  designing  clutter  filters  for 
weather  radar  applications,  as  it  does  not  provide  the  designer  with  sufficient  control 
of  the  filters  stopband  and  transition  region  characteristics. 


(MAGNITUDE) 


STAGGER  RATIO  4:5 
bT  IS  0.26471  ms 
PRF  1  =  944.4  Hz.  PRF  2 


(dB)  vs  DOPPLER  VELOCITY  (m  s) 
FILTER  ID:  24 


10  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3  00000  GHz 
755  5  Hz  MAXIMALLY  FLAT  DESIGN 


Figure  4.6:  |G(-o)  |2  for  MFR  design  -  Ci=4  ,  c2=5  and  N  =10 


4.3  DESIGN  TECHNIQUE:  OUTPUT  SIGNAL  TO  CLUTTER  RATIO 
MAXIMIZATION 


Unlike  the  Coefficient  Selection  and  MFR  design  techniques  which  are  both 
deterministic  approaches  to  the  design  of  Staggered  PRT  clutter  filters,  the  Output 
Signal  to  Clutter  Ratio  Maximization  (OSCRM)  technique  selects  partial  response 
sequences  for  the  Staggered  PRT  in  such  a  manner  as  to  maximize  the  signal  to 
clutter  ratio  (SCR)  in  the  filter’s  output  sequence. 

Given  a  wide-sense  stationary  input  sequence  {x*  >  =  {x  (  dj)},  the  staggered 

PRT  filter’s  output  at  tune  is 

iV-l 

yi  =  £  xi-khij t  .  i  =1  -2 

k  =0 

where  a  +b  denotes  a  modulo  b .  Equivalently,  y/  =  h,HX;  ,  i  =1  - 2  where 
X;  =[  x,  ■  ■  ■  X[ ]T.  h,  =[  ht  0  ■  ■  ■  ht  lV _ i  ]H.  and  superscript  H's  and  T's  indi¬ 
cate  Hermetian  and  standard  transposes,  respectively.  From  equation  3.6  we  see 
that  the  lag  zero  autocorrelation  of  the  filter’s  output  sequence  is 


Ryy  (0)  =  Ex  {  I  y,  I2}  -  |  E.v ,u.„  (h,HX/  (h,HX/ )«}  +  jEXJ  lM  (h,Hx,  ( h,Hx, )»} 


where 


1 
2 

1  ruH 


ho*E x.i  Wc»{x/x;H}h,,  -  h?EXJ  (ldd  (x/  X/H}h, 


R,  - 


^-i-rivXho+h^Rjhj] 


R.xx  Idi  ,0 — 1 Rxx  (dL  .0~di  ,1*  •  ■  Rxx(did)— ^t./V-ll 

Rx*  'i—di  ft)  Rx.x  ^di q  cf,  ] )  .  .  Rxx{di  j-^ , a/ — i ^ 


Rxx^di  jj-i  —di.o)  Rxx  <'dl < d,  q)  .  .  Rxx(di  jj --[—di'N-y) 


(4.4a) 

(4.4b) 


(4.4c) 


Weather  radar  signals  contaminated  by  clutter  may  be  considered  as  the  sum 
of  two  independent  signals  -  the  sampled  returns  due  to  hydrometeors  or  other  parti¬ 
culates  in  the  atmosphere,  {xj^}  ,  and  the  sampled  returns  due  to  clutter  ,  {x£>  .  We 
make  the  assumption  that  both  of  these  signals  may  be  considered  wide-sense  sta¬ 
tionary  over  an  azimuth  dwell  time  (this  assumption  becomes  increasingly  less  valid 
with  increasing  radar  slant  range  due  to  beam  broadening);  and  we  neglect  any 
non-linearities  in  the  radar  receiver  which  would  correlate  the  weather  and  clutter 
signals.  The  output  sequence  of  the  Staggered  PRT  filter  is  then  the  sum  of  the 
independent  contributions  of  the  filtered  weather  and  clutter  processes. 
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Given  the  second  order  statistics  of  the  weather  and  clutter  processes,  the  SCR 
at  ‘he  filter's  output  is  a  scalar  function  of  the  partial  response  sequences  h, .  Denot- 
•ng  the  autocorrelation  matrices  corresponding  to  equation  4.4c  for  the  input  weather 
and  clutter  processes  as  R“  and  R-  ,  respectively,  we  see  that  the  output  SCR  for 
i  =1  +2  is  given  by 


SCR(  h,  )  4 


h,HRfh, 

h.XR'h, 


(4.5) 


Later  on  m  this  section  of  the  report,  we  will  require  models  for  the  weather  and 
clutter  power  spectral  densities  which  are  even  functions  of  Doppler  frequency.  The 
corresponding  matrices  R,a  and  R/  are  then  real,  symmetric  and  positive  semi- 
definite.  The  clutter  power  spectrum  will  be  modeled  as  Gaussian  with  a  non-zero 
spectrum  width.  This  is  sufficient  to  guarantee  that  Rf  will  be  positive-definite  and 
therefore  invertible.  In  order  to  simplify  the  derivation  of  the  OSCRM  technique,  it 
will  be  assumed  for  the  remainder  of  this  section  that  R,a  and  R,°  are  real  and  sym¬ 
metric  and  that  us  invertible.  These  assumptions  in  turn  imply  real  h,  m  the 
context  of  OSCRM. 


The  gradient  of  SCR(  h,  i  is 

v(hTR,-hl  Uh^R/h,  j-vIhT^h,-  )  (  h,T R“ h,  ) 

VSCR<  h,  i  =  — - — - - - — ~2 - — - - 

(hjRrh, 


K  h,  | 

hTR/h, 

l-R/'h,  | 

h,T  R,a  h,  ) 

h‘T 

RjCh, 

2 

Setting  the  preceding  equation  to  zero  and  premultiplying  by  jR,1' j  .  we  obtain 

(ivfVh,  (h^h,  )  =  (h,TR!ahi  )ht- 


or 


Kc\  Rfh, 


h,T  R1a  h. 


>h,-  =  SCR(  h,  )  hj 


(4.6) 


Equation  4.6  is  particularly  interesting  in  that  it  is  an  eigenvalue-eigenvector 
expression  for  the  matrix  |RjC  j  R,a  -  with  eigenvalues  SCR(h,  )  and  eigenvectors 
h,  .  Since  we  arrived  at  this  result  by  evaluating  VSCR(  h,  )  ,  the  desired  partial 
response  sequence  /q  is  the  eigenvector  of  |r,c  j  R,a  corresponding  to  the  eigenvalue 

of  maximum  modulus.  Results  similar  to  this  one,  with  different  derivations,  may 
be  found  in  works  of  Cadzow  [Cadzow80a]  ,  and  Keeping  and  Sure  a  u  [KeepingTSa] 
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The  OSCRM  design  technique  operates  as  follows. 


1)  Select  values  for  Cj  ,  c2  and  the  filter  length  N . 

2)  Select  models  for  the  weather  and  clutter  spectra  which  yield  real 
and  symmetric  R“  and  and  a  positive-definite  R/. 

3)  For  i  =  0,1  .  determine  the  eigenvalues  and  eigenvectors  of  the 
matrix  (r,c  ]  R"  .  The  hL  are  then  the  eigenvectors  corresponding 

to  the  eigenvalues  of  maximum  modulus  for  i  —  0,1  ,  while  the 
clutter  improvement  factors  for  the  assumed  weather  and  clutter 
models  are  the  corresponding  eigenvalues. 

The  use  of  OSCRM  requires  that  some  assumptions  be  made  about  the  statis¬ 
tics  of  the  weather  and  clutter  processes.  For  land-based  weather  radar  there  are 
basically  two  types  of  ground-based  clutter:  scan-modulated  distributed  clutter  such 
as  trees:  and  scan-modulated  discrete  clutter  such  as  isolated  buildings  and  radio 
towers.  At  S-Band,  with  a  one  degree  beam  and  antenna  rotation  rates  on  the  order 
of  3  RPM.  the  spectrum  for  distributed  clutter  is  Gaussian  and  centered  at  zero  fre¬ 
quency  with  a  standard  deviation  of  ~0.25  m/s  ,  while  the  spectrum  for  discrete 
clutter  is  non-Gaussian  and  somewhat  narrower.  The  clutter  model  employed  in  the 
OSCRM  technique  is  the  generally  wider,  distributed  clutter  model.  The  selection  of 
a  model  for  the  weather  spectrum  is  more  difficult.  Although  it  is  Gaussian,  its 
moments  are  a  function  of  the  prevailing  atmospheric  conditions.  For  this  purposes 
of  OSCRM.  we  will  employ  the  generalized  weather  spectrum  depicted  in  Figure  4.7. 

The  inputs  to  OSCRM  are  then  the  stagger  ratio,  N  ,  the  clutter  spectrum 
width,  oc .  and  the  four  breakpoints  of  the  generalized  weather  spectrum, 

iv-\  ,  iv 2  ,  and  w 4.  Of  course.  0  <w  j  <  w2  <  w 3  <  w4  <  —  Hz.  When  OSCRM 

°T 

is  used  to  generate  clutter  niters,  the  values  of  the  parameters  w3  and  w4  are  rela¬ 
tively  unimportant,  as  it  is  the  Doppler  interval  over  which  the  weather  and  clutter 
spectra  coexist  that  drives  the  OSCRM  design.  The  OSCRM  also  includes  a  param¬ 
eter  asc  which  determines  the  SCR  at  the  Staggered  PRT  filter’s  input.  For  a  typi¬ 
cal  clutter  filter  design,  however,  w3—w2  »  oc  and  Wi  »  w2  ~  vc  ,  so  that  asc  has 
fit  tie  effect  on  the  resulting  design. 

Figure  4.8  depicts  an  OSCRM  design  for  a  3  GHz  (S-Band)  radar  with 
<^=0.45455  ms,  corresponding  to  a  Nyquist  velocity  of  55  m/s.  The  stagger  ratio  in 
this  case  is  2:3,  with  a  filter  length  of  N  =3.  Figure  4.8a  displays  the  z-plane  zero 
distribution  for  |G(:c)|2  ,  while  Figure  4.8b  displays  the  corresponding  Doppler 
domain  plot.  The  text  blocks  on  the  lower  right-hand  comer  of  the  plots  display  the 
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Figure  4.7:  Doppler  domain  input  parameters  for  OSCRM  design  technique 
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Z-PLANE  ZERO  PLOT  OF  q  (z>q  (1/z)+C,  (z)C.  (1/z)  FOR  FILTER  233 
NYQUIST  FREQUENCY  IS  (26T  ) Hz 

STAGGER  RATIO  2:3  AUTOCORRELATION-BASED  DESIGN: 

3  FILTER  COEFFICIENTS  INPUT  S/C  POWER  RATIO:  -10  0  dB 

4  ZEROES  AT  z  =  1  OUTPUT  S/C  POWER  RATIOS  62.2  dB.  62.2  dB 

NORMALIZED  CLUTTER  WIDTH:  0  00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS: 

0.00000  0.00200  0  85000  1.00000 

Figure  4.8a:  Z-Plane  zeroes  of  |G(~')|2  for  OSCRM  design  -  Ci=2  ,  c2=3  and 
N= 3 
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(MAGNITUDE)2  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  ID:  233 

STAGGER  RATIO  2:3  AUTOCORRELATION-BASED  DESIGN: 

&T  IS  0  45455  ms  INPUT  S/C  POWER  RATIO  -10.0  dB 

PRF  1  =  1 100  0  Hz,  PRF  2  =  733.3  Hz  OUTPUT  S/C  POWER  RATIOS:  62.2  dB. 

3  FILTER  COEFFICIENTS  62.2  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1  00000 


Figure  4.8b:  |GU’)  |2  for  OSCRM  design  -  Cj— 2  ,  c2=3  and  JV=3 
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OSCRM  parameters  with  ac  and  the  normalized  to  the  Nyquist  velocity.  The 
unnormalized  values  are  crc  =0.15  m/s,  a»1=0  m/s,  m2=0.11  m/s,  w3= 47  m/s  and 
w4= 55  m/s.  These  w were  selected  to  generate  a  generalized  weather  spectrum 
which  essentially  coexists  with  the  clutter  spectrum,  but  with  a  diminished  weighting 
for  the  weather  about  zero  velocity.  osc  ,  the  input  SCR,  was  set  to  -10  dB. 

The  output  SCR’s  for  h0  and  hi  ,  maximum  modulus  eigenvalues  of  the 
(R/pR/"  ,  appear  in  the  lower  right-hand  text  block  as  62.2  dB  for  i  =0,1  ,  indi¬ 
cating  a  clutter  improvement  factor  of  72.2  dB  for  the  assumed  weather  and  clutter 
spectra.  An  equispaced  FIR  filter  with  three  coefficients  would  not  generally  provide 
this  amount  of  clutter  rejection.  Comparing  Figure  4.8a  with  the  corresponding 
MFR  design  in  Figure  3.14a,  we  see  that  the  OSCRM  design  technique  generates  a 
looser  cluster  of  controllable  zeroes  at  z  =1.  This  is  a  result  of  the  OSCRM’s  attempt 
to  cancel  clutter  with  a  non-zero  spectrum  width;  Figure  4.8b  correspondingly 
displays  a  wider  stopband  than  Figure  3.14b.  The  transition  region  and  passband 
characteristics  of  the  design  are  largely  determined  by  the  uncontrollable  zeroes, 
however,  so  that  Figures  4.8b  and  3.14b  agree  in  these  two  respects. 

Figure  4.9  displays  the  clutter  improvement  factor  of  the  design  depicted  in  Fig¬ 
ure  4.8  for  the  distributed  clutter  model  as  a  function  of  the  clutter  spectrum  width 
in  m/ s.  Note  that  the  clutter  rejection  exceeds  40  dB  for  spectrum  widths  up  to  ~0.8 
m/ s.  Also  note  the  minim um  in  the  rejection  characteristic  at  oc  . 

Figure  4.10  displays  the  clutter  improvement  factor  for  scan-modulated  discrete 
clutter  targets  as  a  function  of  antenna  rotation  rate.  Four  data  points  are  plotted. 
The  data  displayed  in  Figure  4.10  were  calculated  with  the  assumption  of  an 
antenna  with  a  uniformly  illuminated  circular  aperture  [Skolnik80a]  .  The  numbers 
next  to  each  datum  in  the  figure  indicate  the  angular  distance  swept  out  to  either 
side  of  the  target  in  the  simulation.  The  filter’s  improvement  factor  for  this  class  of 
clutter  is  acceptable  for  weather  radar  applications. 

Figures  4.11a  and  4.11b  display  the  bias  introduced  by  the  filter  for  R0 
reflectivity  estimates,  as  a  function  of  the  mean  velocity  of  the  weather  spectrum. 
Figure  4.11a  displays  the  bias  for  a  spectrum  width  of  4  m/s,  while  Figure  4.11b 
was  generated  for  a  spectrum  width  of  7  m/s.  Note  that  while  the  performances 
depicted  are  comparable  for  mean  velocities  greater  than  10  m/s,  the  biases  induced 
by  incursions  of  the  weather  spectra  into  the  transition  region  and  stopband  are 
greater  for  narrower  spectrum  widths,  as  expected. 

Figure  4.12  depicts  an  OSCRM  design  which  is  identical  to  that  of  Figure  4.8 
except  that  the  stagger  ratio  is  5:7  rather  than  2:3.  The  point  of  this  figure  is  to 
show  that  the  passband  and  transition  region  characteristics  are  strongly  dependent 
upon  the  choice  of  stagger  ratio.  The  design  displayed  in  Figure  4.12  possesses  a 
transition  region  which  is  half  the  width  of  that  displayed  in  Figure  4.8.  This  is 
achieved,  however,  at  the  expense  of  doubling  the  peak-to-peak  passband  ripple 
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NORMALIZED  TYPE  A  CLUTTER  REJECTION  (ctB)  vs 
CLUTTER  SPECTRUM  WIDTH  (m/s) 

FILTER  ID:  233 

STAGGER  RATIO  2:3  AUTOCORRELATION-BASED  DESIGN: 

ST  IS  0.45455  ms  INPUT  S/C  POWER  RATIO:  -1 0.0  dB 

PRF  1  =  1 100.0  Hz.  PRF  2  =  733.3  Hz  OUTPUT  S/C  POWER  RATIOS.  62.2  dB. 

3  FILTER  COEFFICIENTS  62.2  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0  85000  1.00000 


Figure  4.9:  Distributed  clutter  rejection  for  OSCRM  design  of  figure  4.8 
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NORMALIZED  TYPE  B  CLUTTER  REJECTION  (dB)  vs 
ANTENNA  ROTATION  RATE  (rpm) 

FILTER  ID:  233 


STAGGER  RATIO  2:3 
8t  IS  0.45455  ms 

PRF  1  =  1 100.0  Hz,  PRF  2  *  733.3  Hz 
3  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 


AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10.0  dB 
OUTPUT  S/C  POWER  RATIOS:  62.2  dB, 

62.2  dB 

NORMALIZED  CLUTTER  WIDTH:  0  00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 


ONE-WAY.  HALF-POWER.  ANTENNA  BEAMWIDTH  IS  1 .0° 


Figure  4.10:  Discrete  clutter  rejection  for  OSCRM  design  of  figure  4.8 
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REFLECTIVITY  ESTIMATE  ERROR  (dBz)  vs  MEAN  WEATHER  VELOCITY  (m/s) 

FILTER  ID:  233 


WEATHER  SPECTRUM  WIDTH  IS  4.0  m/s 
MINIMUM  DISPLAYED  VELOCITY:  2.0  m/s 
STAGGER  RATIO  2:3 
8t  IS  0.45455  ms 

PRF  1  =  1100.0  Hz.  PRF  2  =  733.3  Hz 
3  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 
NYQUIST  SPECTRAL  FOLDING  ENFORCED 


Figure  4.11a:  Reflectivity  estimate  bias 
m/s 


AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10.0  dB 
OUTPUT  S/C  POWER  RATIOS:  62.2  dB. 

62.2  dB 

NORMALIZED  CLUTTER  WIDTH:  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 


for  OSCRM  design  of  figure  4.8  -  av 
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REFLECTIVITY  ESTIMATE  ERROR  (dBz)  vS  MEAN  WEATHER  VELOCITY  (m/s) 

FILTER  ID;  233 


WEATHER  SPECTRUM  WIDTH  IS  7.0  m/s 
MINIMUM  DISPLAYED  VELOCITY:  2.0  m/s 
STAGGER  RATIO  2:3 
6T  IS  0.45455  ms 

PRF  1  =  1 100  0  Hz.  PRF  2  =  733.3  Hz 
3  FILTER  COEFFICIENTS 
RADAR  FREQUENCY  IS  3.00000  GHz 
NYQUIST  SPECTRAL  FOLDING  ENFORCED 


Figure  4.11b:  Reflectivity  estimate  bias 
m/ s 


AUTOCORRELATION-BASED  DESIGN: 

INPUT  S/C  POWER  RATIO:  -10.0  dB 
OUTPUT  S/C  POWER  RATIOS:  62.2  dB. 

62.2  dB 

NORMALIZED  CLUTTER  WIDTH:  0.00273 
NORMALIZED  W-SPECTRUM  BREAKPOINTS: 
0.00000  0.00200  0.85000  1.00000 


for  OSCRM  design  of  figure  4.8  -  av=  7 
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(MAGNITUDE)2  (dB)  vs  DOP°LER  VELOCITY  (m/s) 

FILTER  ID:  573 

STAGGER  RATIO  5:7  AUTOCORRELATION-BASED  DESIGN: 

8t  IS  0.45455  ms  INPUT  S/C  POWER  RATIO:  -1 0.0  dB 

PRF  1  .  440.0  H2.  PRF  2  =  314.3  Hz  OUTPUT  S/C  POWER  RATIOS:  46  5  dB. 

3  FILTER  COEFFICIENTS  46  5  dB 

RADAR  FREQUENCY  IS  3.00000  GHz  NORMALIZED  CLUTTER  WIDTH:  0.00273 

NORMALIZED  W-SPECTRUM  BREAKPOINTS 
0.00000  0.00200  0.85000  1.00000 

Figure  4.12:  |G(w)  |2  for  OSCRM  design  -  C]=5  ,  c2=  7  and ^=3 
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relative  to  Figure  4.8.  The  clutter  improvement  factor  for  this  design  is  56.5  dB. 

Within  the  constraints  of  Staggered  PRT  filters,  the  OSCRM  technique  appears 
to  be  the  most  effective  design  technique.  It  allows  the  designer  to  match  the  filter’s 
magnitude  response  to  both  the  weather  and  clutter  spectra;  although  the  resulting 
lag  c  t  and  lag  c2  transfer  functions  are  subject  to  the  phase  difficulties  discussed 
earlier. 


105 


5.0  CONCLUSIONS 

We  have  seen  that  weather  radars  like  the  TDWR  operate  at  S-Band  or  C- 
Band  (2-7  GHz);  and  that  accurate  Doppler  measurements  for  phenomena  such  as 
gustfronts  require  an  unambiguous  velocity  of  at  least  12  m/s.  This  minimum 
unambiguous  velocity  corresponds  to  maximum  unambiguous  ranges  of  approxi¬ 
mately  300  and  150  km  at  S-Band  and  C-Band,  respectively.  We  noted  that  regions 
of  significant  reflectivity  could  be  sensed  beyond  300  km  at  low  elevation  angles,  and 
hence  that  obscuration  of  low  reflectivity,  in-trip,  phenomena  by  high-reflectivity, 
out-of-trip,  weather  would  occasionally  occur.  We  cited  studies  in  which  microburst 
windspeeds  above  30  m/s  had  been  observed,  and  concluded  that  Doppler  weather 
radars  operating  at  S  and  C-Bands  must  contend  with  degrees  of  both  range  obscu¬ 
ration  and  velocity  aliasing.  A  distinction  was  made  between  relative  and  absolute 
velocity  de-aliasing;  it  was  stated  that  absolutely  de-aliased  velocity  observations  are 
required  for  certain  automated  weather  detection  functions  (e.g.  the  Sectorized  Uni¬ 
form  Wind  Algorithm).  We  showed  that  a  clutter  rejection  capability  in  excess  of  50 
dB  is  required  if  a  10  dB  signal  to  clutter  ratio  is  to  be  maintained  for  low- 
reflectivity,  low-altitude,  phenomena. 

We  have  also  seen  that  mean  velocity  estimates  generated  from  a  radar  pulse 
train  with  two  PRT’s  may  be  combined  to  form  an  estimate  which  is  unambiguous 
over  a  velocity  interval  larger  than  that  corresponding  to  either  of  the  two  PRT’s. 
For  two  relatively  prime  positive  integers,  c  l  and  c2,  estimates  from  a  pulse  train 
with  embedded  PRT’s  of  and  c2& r  may  be  combined  to  form  estimates  with  a 

Nyquist  interval  corresponding  to  a  PRT  of  .  Two  such  pulse  trains  were 
described,  the  Batch  PRT  scheme  in  which  the  azimuth  dwell  time  is  divided  into 
two  blocks  of  equispaced  PRT  pulses;  and  the  Staggered  PRT  scheme  in  which  the 
PRT  alternates  from  pulse  to  pulse.  The  Staggered  PRT  scheme  has  the  advantage 
over  the  Batch  PRT  scheme  of  spatial  coherency  for  velocity  estimates  at  the  two 
PRT’s;  and  we  noted  that  the  Batch  PRT  scheme  may  not  effectively  support  velo¬ 
city  de-aliasing  in  the  presence  of  azimuthal  shear. 

5.1  BATCH  PRT 

With  azimuth  dwell  times  of  40  to  50  ms  (20  to  25  degree/ sec  azimuth  scan¬ 
ning  rates),  time-domain  clutter  rejection  techniques  such  as  FIR  filtering  do  not 
seem  practical  with  a  Batch  PRT  sampling  scheme.  This  is  due  to  the  desire  for  a 
narrow  filter  transition  region  and  a  sufficient  number  of  samples  at  the  filter’s  out¬ 
put  to  generate  autocorrelation  estimates  with  small  variances.  A  clutter  rejection 
technique  based  on  the  windowed  periodogram  of  the  radar  return  sequence  was 
described;  and  an  expression  was  developed  which  characterized  the  performance  of 
a  data  taper  in  this  application.  The  taper  performance  was  expressed  as  a  function 
of  the  ratio  of  the  filter  cut-off  frequency  and  the  clutter  spectrum  width.  This  same 
expression  can  be  used  to  characterize  the  performance  of  data  tapers  in  minimizing 
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biases  of  the  spectral  characteristics  of  the  weather  process,  in  particular  the  spec¬ 
trum  width. 

S-Band  clutter  rejection  examples  were  provided  for  the  cases  of  single  PRT  and 
Block  PRT  schemes  with  40  and  50  ms  azimuth  dwell  times.  We  saw  that  the  per¬ 
formance  of  a  given  data  taper  for  this  application  seemed  to  be  driven  by  its 
frequency-domain  sidelobe  characteristics  rather  than  its  main  lobe  width.  For  the 
single  PRT  case,  the  requisite  clutter  rejection  was  available  at  scanning  rates  of  20 
and  25  degrees/ sec  with  a  3  m/s  passband  edge  and  the  Blackman  or  Dolph- 
Chebyshev  tapers.  For  the  Batch  PRT  scheme  with  the  same  range  of  dwell  times 
however,  the  necessary  rejection  was  only  available  with  a  passband  edge  on  the 
order  of  6  m/s,  which  would  not  be  acceptable  for  subsequent  de-aliasing  operations. 
In  the  absence  of  a  substantially  improved  data  taper  design,  we  concluded  that 
either  the  azimuth  dw'ell  time  must  be  increased  beyond  50  ms  or  the  50  dB  rejec¬ 
tion  requirement  relaxed,  if  the  Batch  PRT  scheme  was  to  be  viable. 

5.2  STAGGERED  PRT 

A  filter  structure  capable  of  processing  data  sequences  with  the  unequal  inter- 
sample  delays  embedded  in  a  Staggered  PRT  pulse  train  was  described.  The  filter  is 
comprised  of  two  linear  time-invariant  FIR  filters,  but  the  overall  structure  is  time- 
varying.  Elements  of  a  general  theory  for  these  filters  were  developed,  with  a  partic¬ 
ular  emphasis  on  the  processing  of  stochastic  signals.  One  of  the  distinguishing 
features  of  these  filters  was  that  a  minority  of  the  zeroes  of  the  filter’s  Z-Transform 
were  controllable.  It  was  seen  that  this  filter  structure  possessed  three  effective 
power  transfer  functions,  depending  upon  the  lag  at  which  the  autocorrelation  func¬ 
tion  of  the  filter’s  output  was  computed.  The  power  transfer  function  at  lags 
p  (ci-K2)  •  P  d  was  observed  to  be  real  and  even,  while  the  power  transfer  functions 
at  lags  p  (c  1+c2)+C2  and  p(c1-(-c2)+c1  ,pel  generally  possessed  non-zero  phase 
characteristics.  For  clutter  filtering  applications  with  short  impulse  response  dura¬ 
tions  in  the  component  FIR  filters  ( i.e .  all  controllable  zeroes  near  z  =1),  we  saw  that 
the  selection  of  stagger  ratio  determined  both  the  filter’s  passband  performance  and 
its  z-plane  zero  distribution. 

The  phase  characteristics  of  the  lag  Cj  and  lag  c2  transfer  functions,  and  their 
interaction  with  velocity  estimation,  were  discussed  in  section  3.2.2.  We  saw  that  a 
linear  phase  characteristic  was  possible  for  even  values  of  N,  and  that  the  estima¬ 
tors  of  equations  3.8  could  then  be  employed  to  produce  correctly  de-aliased  velocity 
estimates.  It  was  noted  that  discontinuous  "linear"  phase  characteristics  seem  to  be 
more  common  than  truly  linear  phase  characteristics  and  that  equations  3.8  do  not 
correct  for  phase  discontinuities.  Even  values  of  N  also  resulted  in  discrepancies  in 
the  magnitude  responses  of  the  lag  p(c1+c2)  and  lag  p  (c1-t-c2)+c1  transfer  func¬ 
tions.  We  showed  that  the  magnitude  responses  of  these  transfer  functions  could  be 
matched  perfectly  with  odd  N  and  h  0=h  j. 
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Three  Staggered  PRT  filter  design  techniques  were  presented  along  with 
representative  designs.  The  first  technique,  Coefficient  Selection,  generated  an 
appropriate  eqmspaced  FIR  design  and  converted  that  design  to  a  Staggered  PRT 
partial  response  sequence  by  time-alignment.  We  saw  that  although  the  DFT  is  an 
isometry  -  equivalently,  that  such  a  partial  response  sequence  was  a  minimum 
norm  approximation  to  the  FIR  equispaced  design  in  the  time  and  DFT  domains  - 
the  partial  response  sequences  Z-Transform  will  typically  possess  a  local  maximum 
in  its  relative  error  at  ^=0  for  clutter  filter  designs.  The  second  design  technique, 
MFR,  placed  all  controllable  zeroes  of  the  filter  at  the  origin.  The  resulting  designs 
possessed  a  deep  notch  at  ^-=0  ,  but  MFR  did  not  afford  the  designer  with  adequate 
stopband  and  transition  region  control.  The  final  design  technique  was  the  OSCRM 
technique,  which  maximized  the  filter’s  output  signal  to  clutter  ratio  based  on  spec¬ 
tral  models  for  the  weather  and  clutter  processes.  A  sample  OSCRM  design  was 
presented,  and  its  rejection  performances  for  both  scan-modulated  distributed  and 
discrete  clutter  were  seen  to  be  satisfactory.  Pulse-pair  reflectivity  estimate  errors 
for  this  design  were  also  evaluated. 

The  sample  filter  designs  presented  in  this  report  all  possessed  partial  response 
sequences  with  a  small  number  of  coefficients.  A  number  of  designs  were 
attempted  with  larger  numbers  of  coefficients  (10  to  50)  in  the  hope  that  filters  with 
narrower  transition  regions  could  be  realized.  In  all  these  cases,  the  increased 
passband  ripple  resulting  from  the  large  number  of  uncontrollable  zeros  offsets  any 
performance  benefits  that  might  result  from  a  narrower  transition  region. 

5.3  FURTHER  RESEARCH 

Before  the  Staggered  PRT  scheme  can  be  declared  to  be  either  appropriate  or 
inappropriate  for  Doppler  weather  radar  applications,  an  improved  understanding  of 
the  properties  of  staggered  PRT  filters  and  their  interaction  with  Doppler  moment 
estimation  will  be  required.  The  two  most  outstanding  issues  in  this  respect  are  the 
z-plane  distribution  of  the  filters’  uncontrollable  zeroes,  and  the  phase  responses  of 
the  lag  c !  and  lag  c2  transfer  functions.  In  regard  to  the  former  issue,  the  C,  (z ) 
possess  a  great  deal  of  structure  and  one  would  initially  anticipate  some  success  in 
analytically  characterizing  the  z-plane  locations  as  a  function  of  the  stagger  ratio  - 
especially  for  those  cases  in  which  all  controllable  zeroes  are  placed  at  2=1.  The 
problem,  however,  has  resisted  a  solution.  A  seemingly  promising  approach  would 
be  a  linear  systems  analysis  of  the  staggered  PRT  filter  as  a  single-input,  multi¬ 
output,  system.  There  are  results  available  for  the  eigenvalues  of  sparse  matrices 
which  might  then  be  applied  to  determine  the  zero  locations. 

In  section  3.3.2  we  showed  that  linear  phase  (possibly  with  jump  discontinui¬ 
ties)  was  guaranteed  for  C  0(eJ"^)C\ (eJ~ lT)  with  an  odd  value  of N  and  odd,  in  the 
sense  of  that  section,  fi, .  It  would  be  surprising  if  this  constituted  the  entirety  of  lag 
C}  transfer  functions  with  linear  phase.  If  one  could  select  the  hi  so  that  Cote7"^) 
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and  CpeJ had  the  same  phase  characteristic,  for  instance,  one  would  have  a  lag 
c ,  transfer  function  with  zero  phase.  The  design  techniques  presented  in  this  report 
independently  generated  h 0  and  hx.  A  staggered  PRT  filter  design  technique  which 
simultaneously  generated  the  partial  response  sequences,  incorporating  constraints 
on  the  lag  c  j  and  lag  c2  phase  responses,  would  be  valuable. 

Alternatively,  a  staggered  HR  allpass  filter  could  be  employed  to  adjust  the 
phase  characteristics  of  the  lag  Cj  and  lag  c2  phase  responses.  Preliminary  results 
have  shown  that  such  filters,  applied  in  advance  of  the  partial  response  sequences, 
can  effect,  for  instance,  a  phase  characteristic  with  regions  of  zero  slope.  The  HR 
filters  evaluated  to  date  average  over  input  and  output  samples  separated  by 
(c  1+c2l<h'.  however,  and  thus  possess  Nyquist  intervals  which  are  smaller  than 
those  of  the  overall  Staggered  PRT  filter.  When  C ,)(e~'u^)C  *  (eJu^)  has  linear 
phase,  the  result  is  a  net  lag  cx  phase  response  which  is  periodic  within  the  interval 
(-  1  .  -of1]  ,  with  regions  of  zero  and  non-zero  slope. 

It  was  noted  in  section  3.3.2  that  properly  de-aliased  velocity  estimates  from 
the  staggered  PRT  filter’s  ouput  sequence  could  possess  a  bias.  It  was  suggested 
that  this  results  from  variations  in  the  lag  c}  and  lag  c2  transfer  functions’ 
passband  magnitude  response.  An  alternative  estimation  technique  is  to  use  the 
de-aliased  staggered  PRT  estimate  to  "steer"  an  accurately  aliased  velocity  estimate 
to  the  proper  Nyquist  interval.  This  latter  estimate  would  be  obtained  from  an 
equispaced  clutter  filter  operating  at  a  decimated  sampling  rate,  say  d relative 
to  (V  The  equispaced  filter  could  be  designed  to  have  a  great  deal  less  passband 
variation  than  the  staggered  PRT  filter,  as  well  as  a  narrower  transition  band.  We 
will  refer  to  this  combination  of  equispaced  and  staggered  PRT  filters  as  a  composite 
filter. 

An  equispaced  pulse  sequence  would  be  embedded  in  the  staggered  PRT 
sequence  to  form  the  composite  sequence  of  Figure  5.1,  where  c}=2  ,  c2=3  and 
dc=3.  This  composite  sequence  has  a  period  of  LCM{dc  ,c1+c2)6r  (where  LCM 
denotes  Least  Common  Multiple)  and  a  pulse  density  that  is  slightly  higher  than 
that  of  the  staggered  PRT  sequence  -  in  the  case  of  Figure  5.1,  60%  versus  40%  of 
the  available  pulse  locations  at  PRT  (Sp.  We  note  that  Figure  5.1  includes  some 
pulse  spacings  of  ^  which  would  increase  the  likelihood  of  obscuration.  Figure  5.2  is 
an  overlay  of  an  OSCRM  |G(w)  |2  and  the  squared  magnitude  response  of  a  29 
point  equiripple  clutter  filter  design  (cj=  2  ,  c2=2  ,  N  =3  and  de— 3).  uu  is  55  m/s 
( X=10  cm)  for  this  example  and  the  first  alias  of  the  equispaced  filter’s  stopband  is 
evident  at  2/31^=362/3  m/s.  For  narrow  spectrum  widths  and  mean  velocities  near 
this  aliased  stopband,  only  the  staggered  PRT  filter’s  velocity  estimate  will  be  avail¬ 
able. 

The  velocity  estimation  procedure  is  depicted  in  Figure  5.3.  The  received  pulse 
train  is  decomposed  into  its  staggered  and  equispaced  components.  The  top  chain  in 
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•  COMPOSITE  SAMPLE  SEQUENCE  IS  DECOMPOSED  INTO 
EQUISPACED  AND  STAGGERED  PRT  SEQUENCES 


•  STAGGERED  PULSE  INTERVALS:  C^6T.  C 26T 

•  EQUISPACED  PULSE  INTERVAL:  de6T 


C..  C,.  d.  ARE  INTEGERS 

12  6 


MtI  C26t  | 

-x — x - x — x- 


e.g  .  C,  2.  C,  =  3,  de  -  3 


-X - X- 


-X - X- 


-X - X- 


STAGGERED  PRT 


W. 


Figure  5.1: 


Sample  association  for  composite  pulse  sequence  -  Cj=2  ,  c2=3  and 
<•=3 


no 


’52963  6’ 


COMPOSITE  (MAGNITUDE)2  RESPONSES  (dB)  vs  DOPPLER  VELOCITY  (m/s) 

FILTER  IDs:  STAGGERED  -  233 
EQUISPACED  -  3303 
6t  IS  0.45455  ms 

RADAR  FREQUENCY  IS  3.000  GHz 
STAGGER  RATIOS: 

STAGGERED  FILTER  -  2:3 
EQUISPACED  DECIMATION  RATE  -  3 


Figure  5.2:  Composite  filter  transfer  functions  -  cx=2  ,  c2= 3  ,  de  =3  ,  N  =3  and 

an  equispaced  filter  with  29  coefficients 
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the  filter  is  the  staggered  PRT  velocity  estimation  scheme  described  in  this  report; 
the  bottom  chain  is  conventional  pulse-pair  estimation  at  the  output  of  the 
equispaced  filter.  The  two  velocity  estimates,  labeled  vs  and  uc ,  for  staggered  and 
equispaced,  are  combined  in  the  block  labeled  velocity  estimation  to  produce  the  final 
estimate.  The  mechanism  of  that  block  is  diagrammed  in  Figure  5.4.  The, ratio  of 
the  lag  zero  staggered  and  equispaced  autocorrelation  estimates,  Rs  (0)/i?e(0),  is 
compared  to  a  threshold  K .  If  that  ratio  is  <K ,  the  staggered  PRT  velocity  esti¬ 
mate  is  used  to  determine  the  proper  Nyquist  interval  for  the  equispaced  velocity 
estimate.  If  the  ratio  is  >K ,  the  weather  is  assumed  to  have  aliased  to  the 
equispaced  filter’s  stopband,  and  the  staggered  PRT  velocity  estimate  is  used. 

Figure  5.5  provides  an  example  of  the  velocity  estimate  errors  for  the  composite 
filter  of  Figure  5.2  assuming  infinite  data  and  no  clutter.  |G(^)  |2  was  substituted 
for  the  lag  c ,  and  lag  c2  transfer  functions  -  errors  due  to  the  phase  characteristics 
of  those  transfer  functions  are  not  considered  here.  The  weather  spectrum  width  is  1 
m/s  and  K=  1.1.  The  horizontal  dashed  lines  in  the  figure  identify  velocity  regions 
where  Rs  (0)/Re  (0)  <K ,  the  estimates  in  the  velocity  regions  identified  by  the  hor¬ 
izontal  solid  lines  are  obtained  directly  from  the  staggered  PRT  filter.  As  expected, 
the  staggered  PRT  estimates  are  relied  upon  in  the  region  around  the  alias  of  the 
equispaced  filter’s  stopband.  They  are  also  relied  upon  for  mean  velocities  from  14 
to  32  m/s,  due  to  the  relatively  higher  passband  gain  of  the  staggered  PRT  filter  in 
that  region.  The  magnitudes  of  the  errors  displayed  in  this  figure,  <0.55  m/s,  are 
considerably  less  than  those  of  the  staggered  PRT  velocity  estimates  (e.g.  as  in  Fig¬ 
ure  3.10);  there  are  also  encouraging  preliminary  results  for  the  estimation  of 
reflectivity  and  spectrum  width  with  a  composite  filter.  Composite  filtering  schemes 
merit  further  investigation,  but  improved  control  over  the  phase  responses  of  the 
lags  c  i  and  c2  transfer  functions  will  also  be  necessary  in  this  context. 
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<cc  <oc 


CONDITION: 

Vs  AND  Ve  ARE  BOTH  VALID 

RESULT: 

Vs  USED  TO  DETERMINE  THE  FOLDS 

YES 

IN  Ve.  V  IS  SET  TO  Ve  PLUS  THE 

APPROPRIATE  CORRECTION 

A 

CONDITION: 

Ve  IS  NOT  VALID.  WEATHER  IS 

BEING  ALIASED  TO  A  FILTER  NOTCH 

A  A 

RESULT: 

v  =  vs 

Figure  5.4:  Combining  staggered  and  equispaced  velocity  estimates 
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MEAN  VELOCITY  ESTIMATE  ERROR  (m/S)  vs  MEAN  WEATHER  VELOCITY  (m/s) 

FILTER  IDs:  STAGGERED  -  233.  EQUISPACED  -  3303 

WEATHER  SPECTRUM  WIDTH  IS  1 .0  m/s  ESTIMATOR  IS  BASED  UPON  THE  LAG 
MINIMUM  DISPLAYED  VELOCITY:  3.0  m/s  ESTIMATE  FROM  THE  SELECTED  FILTER  (S). 

6t  IS  0.45455  ms  STAGGERED  TO  EQUISPACED  FILTERS  POWER 

RADAR  FREQUENCY  IS  3.000  GHz  RATIO  SWITCHING  THRESHOLD:  1.10 

STAGGER  RATIOS:  VELOCITY  MODULUS  LIMIT  IS  65.0 

STAGGERED  FILTER  -  2:3  NYQUIST  SPECTRAL  FOLDING  ENFORCED 

EQUISPACED  DECIMATION  RATE  -  3 


Figure  5.5:  Velocity  estimate  errors  for  composite  filter  of  figure  5.2,  aL =1  m/s, 

lag  c  j  and  c2  transfer  functions  replaced  by  |G(u,)  |2 
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Appendix  A:  Gust  Front  Spectrum  Width  Data 

A  gust  front  is  the  leading  edge  of  the  outflow  produced  when  a  thunderstorm 
downdraft  reaches  the  earth’s  surface  and  propagates  horizontally.  Early  and  reli¬ 
able  detection  of  gust  fronts  requires  accurate  Doppler  measurements,  as  a  gust 
front’s  initial  reflectivity  signature  is  often  obscured  by  that  of  its  parent  storm 
[Klingle84a]  .  Table  A.l  displays  seasonal  data  on  gust  front  spectrum  widths  from 
four  regions  East  of  the  Rocky  Mountains.  All  measurements  were  made  within  100 
km  of  the  radar  and  are  in  units  of  m / s.  The  data  displayed  in  the  table  are: 

AvgAvg  - 

per-tilt  average  gust  front  spectrum  widths,  averaged  over  all  tilts 
containing  a  given  gust  front,  then  averaged  over  all  gust  fronts  in 
the  region; 

MaxAvg  - 

per-tilt  maximum  gust  front  spectrum  widths,  averaged  over  all 
tilts  containing  a  given  gust  front,  then  averaged  over  all  gust 
fronts  in  the  region; 

and  Peak  - 

peak  gust  front  spectrum  width  observed  in  a  region. 


Gust  Front  Spectrum  Width  Data 

Region  and 
Source  t 

1 

May-September 

October-Aprii 

MaxAvg 

Peak 

Number  of 
Events 

H 

■ 

N/M 

0 

- 

1 

— 

mm 

7 

n 

5.5 

5 

2.5 

n 

5.5 

NSSL 

[KJingle84a] 

t 

6 

5.6 

■ 

15.0 

4 

6.6 

n.i 

Denver  FL-2 

11 

2.8 

4.0 

6.0 

0 

- 

- 

- 

TABLE  A.l 

t  FL-2  is  the  FAA  S-Band  Doppler  weather  radar  testbed;  NSSL  is  the  National  Severe  Storms 
Laboratory;  N/M  indicates  that  values  were  not  measured  by  the  storm  analyst.  The  FL-2 
data  were  compiled  at  Lincoln  Laboratory  by  Paul  J.  Biron,  Robert  G.  Hallowell,  Mark  A.  Isam- 
inger,  J.  Patrick  Sanford  and  the  author. 
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+  Dr.  Dusan  Zrnic  of  NSSL  indicates  that  these  values  may  have  a  positive  bias  due'to  clutter 
contamination  and  the  fact  that  the  lag  zero  autocorrelation  estimates  used  in  computing 
NSSL’s  spectrum  width  display  are  not  col  lected  for  system  noise. 

The  Memphis,  Huntsville  and  Denver  data  are  consistent  in  all  three  categories 
with  an  average  AvgAvg  value  of  2.5  m/s  and  average  MaxAvg  value  of  4.1  m/s. 
The  overall  peak  value  for  those  regions  is  6  m/s.  The  magnitude  of  the  bias,  if 
any,  in  the  NSSL  data  is  unknown.  And,  if  the  NSSL  measurements  are  accurate, 
such  energetic  gust  fronts  would  likely  contain  high  amounts  of  debris  and  hence 
possess  prominent  reflectivity  signatures.  We  propose  that  reflectivity  measure¬ 
ments  along  with  accurate  Doppler  measurements  for  weather  of  spectrum  widths 
up  to  6  m/s  will  be  sufficient  for  the  detection  of  gust  fronts. 
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